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HIGGS PARTICLES 
Introduction and Summary 

A. Djouadi^'^ and P.M. Zerwas^ 

^ Deutsches Elektronen-Synchrotron DESY, D-22603 Hamburg, FRG. 
^ Institut fxir Theoretische Physik, Universitat Karlsruhe, D- 76128 Karlsruhe, FRG. 

The search for scalar Higgs particles and the exploration of the mechanism which 
breaks the electroweak symmetry, will be one of the major tasks at future high-energy 
e'^e~ colliders. In previous studies it has been shown that e^e~ linear colliders operating 
in the energy range \/i ~ 300 to 500 GeV with a luminosity of / £ ~ 20 fb~^ are ideal 
machines to search for light Higgs particles. 

In the Standard Model (SM) the Higgs mass range Mh ^ 200 GeV is easy to cover at 
these energies. This intermediate Higgs mass range is one of the theoretically most favored 
ranges, allowing the particles to remain weakly interacting up to the GUT scale A ~ 10^^ 
GeV [a prerequisite for the perturbative renormalization of the electroweak mixing angle 
from the GUT symmetry value 3/8 down to the experimentally observed value at low 
energies]. The search of intermediate- mass Higgs bosons can be carried out in three 
different channels: the Higgs-strahlung process e+e" — » ZH and the fusion mechanisms 
WW/ZZ — > H . The cross sections are large and the properties of the Higgs boson, in 
particular spin-parity quantum numbers and couplings to gauge bosons and fermions, can 
be thoroughly investigated, allowing for crucial tests of the Higgs mechanism. 

In the Minimal Supersymmetric Standard Model (MSSM) the Higgs sector is extended 
to three neutral scalar and pseudoscalar particles h/H, A and a pair of charged particles 
H"^. The lightest Higgs boson h has a mass Mh < 140 GeV and can be detected in the 
entire MSSM parameter space either in the Higgs-strahlung process, e+e" hZ, or by 
the complementary mechanism of associated production with the pseudoscalar particle, 
e'^e~ — > hA. Moreover, there is a substantial area in the MSSM parameter space where 
the heavy Higgs bosons can be also found; for a total energy of 500 GeV this is possible 
if the H, A and masses are less than 230 GeV. Similar to the SM, various properties 
of these Higgs bosons can be investigated. 

Higher energies are required to sweep the entire mass range of the SM Higgs particle, 
Mfj < 1 TeV. The high energies will also be needed to produce and to study the heavy 
scalar particles in extensions of the SM, such as the MSSM, if their masses are larger 
than ~ 250 GeV. Masses of the heavy Higgs bosons in this range are suggested by grand 
unified supersymmetric theories. In e"'"e~ collisions, these experiments can be performed 
in the second phase of the colliders with a cm. energy up to 1.5 to 2 TeV. In this report, 
we analyze the potential of a 1.5 TeV e+e^ linear collider, with an integrated luminosity 
of / £ ~ 200 fb~^ per annum to compensate for the drop of the cross sections at high 
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energies. We will study the heavy Higgs particles in the Standard Model, in the minimal 
supersymmetric extension and in other more speculative scenarios. 

In the Standard Model, the main production mechanisms of Higgs particles, e^e~ 
HZ and e+e^ — > uuH/e^e^H, will be discussed and the cross sections, including the 
interference between the Higgs-strahlung and the fusion processes, will be given [2] . The 
double Higgs production process, in which the trilinear Higgs coupling can be determined 
[therefore leading to the first non-trivial test of the Higgs potential], will be discussed 
in Ref. [3]. Finally, possible effects of New Physics beyond the SM on production cross 
sections and angular distributions of Higgs bosons, will be summarized in Ref. [4]. Conse- 
quences of a model in which the Higgs boson interacts strongly with scalar singlet fields 
in a hidden sector, are described in Ref. [5]. 

Subsequently, we will investigate the properties of the heavy Higgs particles in super- 
symmetric extensions of the SM. We will restrict ourselves first to the minimal extension 
which is highly constrained, parameterized by only two free parameters at the tree-level: 
a Higgs mass parameter [generally the mass the pseudoscalar Higgs boson M^] and the 
ratio of the vacuum expectation values of the two doublet fields responsible for the symme- 
try breaking, tg/3, which in grand unified supersymmetric models with Yukawa coupling 
unification is forced to be either small, tg/3 ~ 1.5, or large, tg/3 ~ 50. 

The various decay modes of the heavy CP-even Higgs boson H, the pseudoscalar 
boson A and the charged Higgs particles will be analyzed in Ref. [6] , in particular the 
decays into supersymmetric particles, charginos, ncutralinos, squarks and sleptons. The 
production of the heavy Higgs particles, primarily through the channels e^e~ HA and 
H^H~ , will be also discussed in Ref. [6] and the complete one-loop electroweak radiative 
corrections of the cross sections will be summarized [7,8]. We will finally discuss the 
multiple production of the SM and the fight MSSM Higgs bosons in Refs.[9,10]. Some of 
these processes will allow us to determine the fundamental Higgs trilinear couplings. 

A brief discussion of the Higgs sector in the next-to-minimal supersymmetric exten- 
sion of the Standard Model, Ref. [11], concludes this report. 
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Higgs-strahlung and Vector Boson Fusion 
in e+e Collisions 

W. KiLiAN, M. Kramer, and P.M. Zerwas 
Deutsches Elektronen-Synchrotron DESY, D-22603 Hamburg/FRG 

Abstract 

Higgs-strahlung e~^e~ ZH and WW {ZZ) fusion e+e^ i^e^eH {e^e^H) are the 
most important mechanisms for the production of Higgs bosons at future e~^e~ hnear 
cohiders. We have calculated the cross sections and energy /angular distributions of 
the Higgs boson for these production mechanisms. When the Z boson decays into 
(electron-)neutrinos or e~^e~ , the two production amplitudes interfere. In the cross- 
over region between the two mechanisms the interference term is positive (negative) 
for z^efg {e~^e~) decays, respectively, thus enhancing (reducing) the production rate. 



The analysis of the mechanism which breaks the electroweak gauge symmetry SU (2)l x 
U{1)y down to U{1)em, is one of the key problems in particle physics. If the gauge 
fields involved remain weakly interacting up to high energies - a prerequisite for the 
(perturbative) renormalization of sin^ 9w from the symmetry value 3/8 of grand- unified 
theories down to a value near 0.2 at low energies - fundamental scalar Higgs bosons 
must exist which damp the rise of the scattering amplitudes of massive gauge particles 
at high energies. In the Standard Model (SM) an isoscalar doublet field is introduced 
to accomodate the electroweak data, leading to the prediction of a single Higgs boson. 
Supersymmetric extensions of the Standard Model expand the scalar sector to a spectrum 
of Higgs particles 0. The Higgs particles have been searched for, unsuccessfully so far, 
at LEPl, setting a lower limit on the SM Higgs mass of mn > 65.2 GeV 0. The search 
for these particles and, if found, the exploration of their profile, will continue at LEP2 [^, 
the LHC 0, and future e~^e~ linear colliders 0. 



Figure 1: Higgs-strahlung and vector boson fusion of (CP-even) Higgs bosons in e~^e 
collisions. 

In this note (see also 0) we will focus on the production of scalar Higgs bosons in e+e" 
collisions. The main production mechanisms for these particles are Higgs-strahlung 
and WW {ZZ) fusion iPHTll] [supplemented in supersymmetric theories by associated 
scalar /pseudoscalar Higgs pair production]. In particular, we will present a comprehensive 
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analysis of the interplay between the production mechanismsj^ (Fig.|ip 



Higgs-strahlung 
WW fusion 
ZZ fusion 



e+e" ^ ZH ^ uvH {e+e~H) 
e~^e~ — i> e~^e~H 



(1) 



For z/gZ/g and e'^e~ decays of the Z bosons, the two production amplitudes interfere. 
It turns out that the interference term is positive for i/gi^e and negative for e^e~ decays, 
respectively, in the cross-over region between the two mechanisms. The interference effect 
had been noticed earlier PJTB[|; however, we improve on these calculations by deriving 
analytic results for the energy and polar angular distribution of the Higgs particle {Eh, 0) 
in the final states of e^e~ — > if + neutrinos and e^e~ He^e~ . This representation can 
comfortably serve as input for Monte Carlo generators like PYTHIA/JETSET and 
HZHA which include the leading QED bremsstrahlung corrections and the important 
background processes. 



Total cross sections. The cross section for the Higgs-strahlung process can be written 
in the following compact form: 

a e+e- ^ ZH) = (vl + al) ^f- 2 

where -y/i is the center-of-mass energy, and = —1, = —1+ 4sin^6'vi/ are the Z 
charges of the electron; A = (1 — {tjih + rnzY/s) (1 — {niH — nizY/s) is the usual two- 
particle phase space function. So long as the non-zero width of the Z bosonQ is not taken 
into account, the cross section rises steeply at threshold ~ (s — {ttih + mzYY^"^- After 
reaching a maximum [about 25 GeV above threshold in the LEP2 energy range] , the cross 
section falls off at high energies, according to the scaling law ~ Qwl^ asymptotically. 
Thus, Higgs-strahlung is the dominant production process for moderate values of the 
energy. The cross section (||) for Higgs-strahlung is reduced by a factor 3 x BR,^ = 20% 
if the final state of Z decays is restricted to neutrino pairs. 

The total cross section for the WW {ZZ) fusion of Higgs particles can be cast into a 
similarly compact formQ 



a 



Glm\, /-i ^ /-i dy 



dx 



QAV^tt^ JxH Jx [l + {y-x)/xvf 



v' + a'yf{x,y) + Av'a'g{x,y) (3) 



^ We will concentrate first on the Standard Model (SM) ; the extension to the Minimal Supersymmetric 
Standard Model (MSSM) is trivial as will be demonstrated in the last section of this note. 



The results presented in this note are insensitive to non-zero width effects of the Higgs boson |16 
For SM Higgs masses below 100 GeV, Th is at least three orders of magnitude smaller than Tz; for larger 
Higgs masses, mn can be reinterpreted as the effective invariant mass of the Higgs decay products. 

•^The variable x is the invariant mass squared of Ve plus H , {x — y) the 4-momentum transfer squared 
from to (all momenta in units of the total energy). 
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where V denotes either W or Z, the charges are v = a = \pl {v = Ve and a = a^) for 
WW (ZZ) fusion, respectively, and 




log(l + z] 



+ 



X z^{l — y) 
y'^ 1 + z 



with xh = rnji/s, xy = rriy/s and z = y{x — xh) / {xxy)- For moderate Higgs masses and 
energies, the cross section, being 0{g^r), is suppressed with respect to Higgs- strahlung 
by the additional electroweak coupling. The smaller value of the Z-electron coupling 
suppresses the ZZ fusion process by an additional order of magnitude compared to WW 
fusion. At high energies, the WW fusion process becomes leading, nevertheless, since the 
size of the cross section is determined by the W mass, in contrast to the scale-invariant 
Higgs-strahlung process. 



cr e e 



w 



1 + 



m 



H 



log- 



2 1 



m 



H 



m 



H 



G'pTnfj/ s 
log- 



Ay/2 



m 



(4) 



H 



The cross section rises logarithmically at high energies, as to be anticipated for this t- 
channel exchange process. 



Differential cross section and interference for WW fusion. The compact form (|^) 
for the fusion cross section cannot be maintained once the interference term between vector 
boson fusion and Higgs-strahlung is included. Moreover, since in the case of WW fusion 
the integration variables x and y do not correspond to observable quantities, the formula is 
useful only for calculating the total cross section without experimental cuts. Nevertheless, 
similarly compact expressions can be derived in this general case by choosing the energy 
Eh and the polar angle 9 of the Higgs particle as the basic variables in the e^e~ cm. 
frame. The overall cross section that will be observed experimentally for the process 



e e 



H 



receives contributions 3 x ^5 from Higgs-strahlung with Z decays into three types of 
neutrinos, from WW fusion, and Qi from the interference term between fusion and 
Higgs-strahlung for VgVf. final states. We findQ for energies \/i above the Z resonance: 



do(Hvv) G\m\'p 



dEf{ d cos 9 



V2 7l^S 



(Sgs + Gi + Gw) 



(5) 



*The analytic result for Qw had first been obtained in Ref . |1 1 
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Figure 2: Total cross sections for the processes e~^e~ Huu and e+e" — > He^e~ as 
a function of the Higgs mass. The cross sections are broken down to the three com,po- 
nents Higgs-strahlung, vector boson fusion, and the interference term, "thr" denotes the 
maximum Higgs mass for on-shell ZH production, "tot" is the total cross section. In 
e^e" — > HVi/ (above) the interference term is negative for small Higgs masses, for large 
Higgs masses positive. In e+e" He^e~ (below), the interference term is of opposite 
sign. 
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with 



Qs 
Qi 



96 



SSp + Si,S2 



[s - mlf [{s^ - m|)2 + m|r|] 



(6) 



(f e + Oe) cos'' 6 



W 



{s - m|) [{s^ - m|)2 + m|r|] 

1 



2 - (/ii + 1) log - (/i2 + 1) log 



X 

cos^ 9w 
siS2r 



hi - 1 

(/ll + l)(/l2 + l) 



/i9 



+ (/ll + l)(/l2 + l) 



(7) 



2ti 



+ 



2t, 



h2-l hi -I 



hi -I 



+ 



hi 



t2 + S' 



6s 



/3tit2 

V r 



(8) 



The cross section is written explicitly in terms of the Higgs momentum p = \J Efj — m|^, 
and the energy = y/s — Eh and invariant mass squared Sy = el — of the neutrino 



pair. In addition, the following abbreviations have been adopted from Ref. |11 



■51,2 = \^{^v ±pCOs6) 

/ii,2 = 1 + 2m^/si,2 

Cx = 1 - 2sSu/{SiS2) 



X X 

To derive the total cross section a{e~^e~ 
be integrated over the region 



tl,2 — ^1,2 + Cx^2,l 

r = h\ + h\ + 2cy,hih2 — 

hih2 + + 

^ = log , , , — r 

hih2 + c^ - 

—>■ Huu), the differential cross section must 



— 1 < cos 6' < 1 and uih < Eh < 



1 + 



m 



H 



(9) 



Differential cross section and interference for ZZ fusion. Similarly, the overall 
cross section for the process 

e~^e~ H + e~^e~ 

receives contributions Qs from Higgs-strahlung with Z decays into electron-positron pairs, 
Qz from ZZ fusion, and Qj from the interference term between fusion and Higgs-strahlung: 



with 
Gs 

01 



da{He~^e-) _ G%m%p 
dEnd cos 9 ^/2n^s 



{Gs + Gi + Gzi + Gz2) 



(10) 



SSe + S1S2 

192 (s-m|)'[(se-m|)2 + m|r|] 
= + alf + Avlal 



Se - mi 



64 



[s - ml) [(se - m|)2 + m|r|] 



2-{h, + 1) log - ih + 1) log + ih + m, + 1)4 



32 51^2^ 

2ti 



^(/ll + l)(/l2 + l) 
2t2 



2 2 

+ 



64 



/if-1 /li 



3tit2 



(12) 



C 



h2-l hi-1 



+ (ti + t2 + 4') — 



t/Z2 = -TT (1 - Cx 

16siS2r 



2 2 

+ 



6s 



h\-\ hl-1 



3tlt2 



C 



(13) 
(14) 



where the same abbreviations as in the formulae following Eq.(§), with the appropriate 
replacements u —>■ e and W ^ Z, have been used. 



To interpret the results, we display the three components of the total cross sections 
a{e'^e~ —>■ Hvv) and (j{e^e~ —>■ He^e~) in Fig.^ for the linear collider energy y/s = 500 
GeV in the cross-over region.^ 

While the energy distribution of the Higgs particle peaks at Eh ~ (s + m|^ — m|)/2i/s 
for Higgs-strahlung, it is nearly flat for WW fusion (Fig.^, left). Only with rising total 
energy the lower part of the Higgs spectrum becomes more pronounced. The angular dis- 
tribution for Higgs-strahlung is almost isotropic at threshold while the standard sin^ 9 law 
is approached, in accordance with the equivalence principle, at asymptotic energies (Fig.^, 
right). The angular distribution peaks, by contrast, in the WW fusion process at — 
and TT for high energies as expected for t-channel exchange processes. 



Polarized beams. At linear colliders the incoming electron and positron beams can be 
polarized longitudinally. Higgs-strahlung and WW fusion both require opposite helicities 
of the electrons and positrons. If cu^l^^rl denote the cross sections in e'^e~ — >■ Huu for 
unpolarized electrons/positrons, left-handed electrons/right-handed positrons, and right- 
handed electrons/left-handed positrons, respectively, we can easily derive, in the notation 
of Eq.(l): 

au oc 3Qs + Qi + Gw (15) 
(Tlr cx QOs + ^Gi + ^Gw (16) 
(JRL OC 6Qs (17) 

The cross section for WW fusion of Higgs particles increases by a factor four, compared 
with unpolarized beams, if left-handed electrons and right-handed positrons are used. By 
using right-handed electrons, the WW fusion mechanism is switched off. [The interference 
term cannot be separated from the WW fusion cross section.] 

^Note that Higgs-strahlung dominates WW fusion at 500 GeV for moderate Higgs masses only if the 
total ZH cross section is considered. 
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For the process e~^e He^e , the pattern is shghtly more comphcated: 

CTu oc gs + Gi + Qzi + Qz2 (18) 

oc ^yr^Qs + 2 ^ ^ (^, + ^.0 (19) 

(t;2 + a2) [vl + a2)2 + 4t;2a2 

O-LL = O-RR (X 2Qz2- (21) 

However, since <^ a^, the difference between aRL and 0"^/? is suppressed. 

Supersymmetric CP-even Higgs bosons. It is trivial to transfer all these results 
from the Standard Model to the Minimal Supersymmetric Standard Model (MSSM). Since 
the couplings to Z/W gauge bosons in the MSSM are shared by the CP-even light 



and heavy scalar Higgs bosons, h and H, respectively, only the overall normalization of 
the cross sections is modified with respect to the Standard Model: 

c^(/i)MSSM = sin^(/5-a) X a(i/)sM (22) 
cr{H)MSSM = cos^{p-a)xa{H)sM (23) 

Higgs-strahlung, vector boson fusion, and the interference term are affected in the same 
way. [The angle a is the mixing angle in the CP-even Higgs sector while the mixing angle 
f3 is determined by the ratio of the vacuum expectation values of the two neutral Higgs 
fields in the MSSM. A recent discussion of the size of the coefficients sin^ / cos^(/3 — a) 



may be found in Ref.[21 . 
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Figure 3: Energy distribution (left) and angular distribution (right) of the Higgs bosons 
for the three components of the cross section /Hs = Higgs-strahlung; WW = fusion; intf 
— interference term]. The individual curves are normalized to the total cross section. The 
Hs peak extends up to maximal values of 0.22 GeV""*^. The total cross section is 69.4 fb. 
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Associated Pair Production of the SA4 Higgs 
and the Probing of the Higgs Self-CoupUng 
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E-mail: ctLopin@lapplLpO.in2p3.fr 



1 Associated Higgs Pair Production 

The interest in double Higgs production is the probing of the triple Higgs self coupling. It 
has been considered in e~^e~ sometime ago|T|. The most efficient means for double Higgs 
production is e^e^ —>■ Uei^eHH (see fig. Q). Double Higgs bremstrahlung (e"^e^ —>■ ZHH) 
is only competing at relatively low energies where the event sample is too low to be useful. 
The equivalent loop-induced double Higgs production in e+e" has been found to be much 
too small 1^ and is not sensitive to the coupling. However, the 77 mode can form a 
= state and therefore 77 HH is a candidate for testing the coupling [Q]. It 
has been pointed out recently that another interesting process is 77 W^W~ HH p| 
that is expected to compete with double Higgs production in e"'"e~ . The reason is that 
in the TeV range, W fusion processes are very much enhanced. The sub-process involved 
is W^W~ —>■ HH , where the dominant helicity amplitude is: 



^ [ PhPw \x - Xq X + Xq 



H , 02 , ^4 



Where h^, is the anomaly in the triple Higgs coupling g, i.e. g = h^gsm where gsm is 
the minimal standard model coupling of H^. We also denote (3w,h = y^l — AM^^^j^/ s, 
r = Ml/M^, xo = (1 + /3],)/2f3w/3H, x = cosO. 

Figure ^ shows that at 2TeV, the cross sections drops precipitously with increasing 
Higgs mass. One can also notice that in 77 —>■ W^W~ HH the external outgoing W 
(mainly transverse) are produced at small angle and take a large amount of energy. For 
the fusion diagrams of this process, the internal W triggers W^W~ — > HH , which 
implies that these diagrams dominate for a heavy Higgs. When convoluting with the 
much advertised photon spectra of[|| and for = lOOGeV^ the cross section of 77 
W^W~ HH drops by about at least a factor of 2 compared with the result without 
convolution. For W^-fusion-like processes, the internal W^s are almost on shell and one 
may wonder if some structure functions could reproduce the exact results. 
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Figure 1: Comparison of cross sections for double Higgs production at e+e and 77 reactions 
for a light Higgs Mh = lOOGeV. 



Figure 2: Higgs mass dependence of the 77 W~^W~ HH cross section at 2TeV. The contri- 
bution of the diagrams involving the triple Higgs vertex (Signal) and the rest (Background) is 
shown separetely. Note the strong interference that occurs in the SAi especially for large Mh. 



15 



Figure 3: Comparing the result of the Wl effective approximation [a ) to the exact result 
^exact £qj, g+g- i/ei'eHH (left) and 77 — > W^W^ HH (right) for a light Higgs and a heavy 
Higgs. Also shown is the asymptotic analytical cross section a^^^. a^'^ is the cross section 
with both outgoing VF's transverse. 



2 The structure function approach 

There have been numerous derivations of the distribution (or structure function) of the W 
inside the light fermions (quarks and electrons) 0. For the effective W approximation, 
the most interesting aspect concerns the Wl content, which has been used to investigate 
manifestations of models of symmetry breaking and Higgs production. The Wl distribu- 
tion inside the photon has only very recently been studied [0. For the case of the heavy 
Higgs the approximation is excellent, already at 2TeV. However, for a light Higgs, the 
approximation is not good and reproduce only the energy behaviour. If one makes the 
further approximation that the hard process cross section is independent of the energy, 
this additional "asymptotic" approximation only reproduces the energy behaviour as well 
as the order of magnitude (even for a heavy Higgs, see figs. ^. 
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Figure 4: The distribution in the reconstructed angle 9* for the signal, background and the 
interference in the case of 77 W~^W~HH without convolution with photon spectra. 

3 Identifying and measuring the Higgs triple vertex 

There is a specific signature of the coupling in all processes that we have studied. Once 
we note that the two Higgses that originate from this vertex can be regarded as produced 
by a scalar H* then in the centre of mass system of the pair, the angular distribution of 
the Higgses is flat. Therefore, we suggest to reconstruct the angle, 9*, measured in the 
centre-of-mass of the pair, between the Higgs direction and the boost axis or the direction 
of the beam. For the "signal"]^ the distribution is flat, while the "background" is peaked 
in the forward/backward direction (see fig. ^). We therefore consider the ratio R of events 
that verify | cos(^*)| < cos(6'q)^ over the number of events outside this region. Assuming 
a total integrated luminosity of 300/6~^, and a 50% efficiency for the reconstruction of 
the double Higgs events one obtain ~ 68 e+e^ Ve^eHH events for Mh = lOOGe^. 
Here the criterion for detection of an anomaly in is a 50% deviation in the expected 
number of events, provided one has at least 30 events. We conclude that with the total 

^We call "signal" the part of the amplitude that include the vertex. The rest is called "background" . 
^ COS0Q = 0.5 is taken in the following. 
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Figure 5: Dependence of the ratio R on h^. 



cross section one would only be able to claim New Physics if dh^ < —0.75 or dh^ > 2. 
For Mh = 4006*61^, SAi values will not lead to a measurement, however if {Sh^l > 1 a 
signal will be recorded (with more than 30 events) and would be a clear indication for an 
anomalous coupling. For = lOOGeV^ where one has enough events for a SAi value, 
the ratio R is much more powerful in constraining the coupling. First the event sample 
within I cos^*| < 0.5 is about 7 out of 60 outside this region. Assuming that the ratio can 
be measured at 20%, we find —.10 < {Sh^l < .15 (see fig. H) which means a precision of 
about 10% on h^. For 77 — > W~^W~ HH and considering the effective 77 luminosity, for 
Mh = lOOGeV, one can hope to collect 15 events. In view of this number the criterion 
for detection of non-standard values is 100% deviation in the number of events. However, 
for a Higgs mass of 400GeV the effect of an anomalous coupling are dramatic and, by 
far, much more interesting than in e~^e~ . Requiring observation of at least 15 events for 
Mh = 400GeV (within the SA4 one expects only 3) useful constraints on the coupling 
can be set: —.7 < Sh-^ < 0.5. There is thus a complementarity between the e^e~ and 
the 77 depending on the Higgs mass in probing the Higgs triple vertex. As for the ratio 
R, taking Mh = lOOGeV it is unlikely that with the number of total WWHH events 
at 77 one would be able to make such a measurement, nonetheless even if this ratio 
were measured with the same precision as in e^e~ one would not constrain the couplings 
further than what is achieved in the classic e^e~ mode. For Mh > QOOGeV, 77 HH 
is the only reaction where useful limits can be set. Thus, there is at a 2TeV collider a 
very nice coverage of the sensitivity by all three reactions. 



18 



4 Conclusions 



We have seen that a 2TeV e+e^ colhder with the reahstic luminosities expected for this 
machine one may hope to achieve a measurement of the tri-hnear couphngs at the level of 
10% (for a light Higgs). The results are also encouraging in the sense that the e^e~ and 
the 77 modes can cover different ranges of the Higgs mass. We find that for a light 
Higgs (up to 250GeV) the best limits on the couplings come from e^e~ —>■ v^^eHH . 
However, for heavier Higgses up to mass of SOOGey, the best channel is the associated 
double Higgs production in 77 . For still heavier masses, the one-loop induced 77 HH 
is by far better. The variable R clearly helps in discriminating the triple Higgs vertex. 
As a by-product we have verified the validity of the distribution function describing the 
longitudinal W content of the photon. We should also insist on the complementarity of 
the e^e~ and the 77 modes of the linear collider for these studies. 
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Anomalous Couplings in the Higgs-strahlung Process 

W. KiLiAN, M. Kramer, and P.M. Zerwas 
Deutsches Elektronen-Synchrotron DESY, D-22603 Hamburg/FRG 

Abstract 

The angular distributions in the Higgs-strahlung process e"'"e~ HZ Hff are 
uniquely determined in the Standard Model. We study how these predictions are 
modified if non-standard couplings are present in the ZZH vertex, as well as lepton- 
boson contact terms. We restrict ourselves to the set of operators which are singlets 
under standard SU3XSU2XU1 transformations, CP conserving, dimension 6, helicity 
conserving, and custodial SU2 conserving. 

The Higgs-strahlung process [l|] 

e+e- ^HZ^ Hff (1) 

together with the WW fusion process, are the most important mechanisms for the produc- 
tion of Higgs bosons in e+e~ collisions Since the ZZH vertex is uniquely determined 
in the Standard Model (SM), the production cross section of the Higgs-strahlung process, 
the angular distribution of the HZ final state as well as the fermion distribution in the Z 
decays can be predicted if the mass of the Higgs boson is fixed 0]. These predictions may 
be modified when deviations from the pointlike coupling are present, which can occur in 
models with non-pointlike character of the Higgs boson itself or through interactions be- 
yond the SM at high energy scales. Since the effective energy scale of the Higgs-strahlung 
process is set by the cm. energy y/s, while the fusion processes are essentially low-energy 
processes with an effective scale of the order M^, new interactions manifest themselves 
more clearly in the total cross section and angular distributions for the Higgs-strahlung 
process (see also 0). 

Operator basis. Deviations from the pointlike coupling can occur in models with non- 
pointlike character of the Higgs boson itself or through interactions beyond the SM at 
high energy scales. We need not specify the underlying theory but instead we will adopt 
the usual assumption that these effects can globally be parameterized by introducing a 
set of dimension-6 operators 

^ = ^SM + E ^(^^ (2) 

i 

The coefficients are in general expected to be of the order 1/A^, where A denotes the energy 
scale of the new interactions. However, if the underlying theory is weakly interacting, the 
ai can be significantly smaller than unity, in particular for loop-induced operators. [It is 
assumed a priori that the ratio of the available cm. energy to A is small enough for the 
expansion in powers of 1/A to be meaningful.] 
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If we restrict ourselves to operators [^] which are singlets under SU3 x SU2 x Ui 
transformations of the SM gauge group, CP conserving, and conserving the custodial SU2 
symmetry, the following bosonic operators are relevant for the Higgs-strahlung process: 

Oa, = ^|9,(^V)p (3) 
O^w = \v^Wl^ (4) 
O^B = \v^Bl^ (5) 

where the gauge fields W^,B are given by the Z, 7 fields. This set of operators is par- 
ticularly interesting because it does not affect, at tree level, observables which do not 
involve the Higgs particle explicitly. [It is understood that the fields and parameters 
are (re-)normalized in the Lagrangian C in such a way that the particle masses and the 
electromagnetic coupling retain their physical values.] 

In addition, we consider the following helicity-conserving fermionic operators which 
induce contact terms contributing to e^e~ —>■ ZH: 

Oli = {^HD^^){hriL) + h.c. (6) 
Ol3 = (<^V«zD^^)(4rV^L)+h.c. (7) 
Or = {y^UD^y^){eRj^eR)+h.c. (8) 

[£l and cr denote the left-handed lepton doublet and the right-handed singlet, respec- 
tively. The vacuum expectation value of the Higgs field is given by {(p) = {0,v/^/2) 
with V = 246 GeV, and the covariant derivative acts on the Higgs doublet as D^^ = 
— ^gT°'WI^ + ^g'Bfj^.] Helicity- violating fermionic operators do not interfere with the 
SM amplitude, so that their contribution to the cross section is suppressed by another 
power of A^. The helicity-conserving fermionic operators modify the SM Zee couplings 
and are therefore constrained by the measurements at LEPl; however, it is possible to 
improve on the existing limits by measuring the Higgs-strahlung process at a high-energy 
e~^e~ collider since the impact on this process increases with energy |^. 



Figure 1: Anomalous ZZH/'yZH couplings and e^e ZH contact terms in the Higgs- 
strahlung process. 

The effective ZZH and the induced 'jZH interactions (Fig.|l], left diagram) may be 
written 

CzzH = gzMz{^-^Z,Z^H+'^Z,,Z^^H^ (9) 
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C^ZH = gzMz^-^Z^^A^^H (10) 



where gz = Mz \JA^/2Gf. Additional operators Z^,Z^''^yH and Z^A^'''^yH are redundant 
in this basis: They may be ehminated in favor of the other operators and the contact terms 
by applying the equations of motion. The remaining coefficients are given by 

ao = -]^aQ^v'^/J\^ (11) 

ai = 4:gz'^ {c^a^pw + s^^a^ps) / J^^ (12) 

hi = Agz'^cwSwi-Oi^w + oi^b) /^^ (13) 

where sw and cw are the sine and cosine of the weak mixing angle, respectively. 

In the same way the eeHZ contact interactions (Fig.|l], right diagram) can be defined 
for left /right-handed electrons and right /left-handed positrons 

^eeZH = gzMz [cLCL^eLH + CRCR^eRH] (14) 

with 

cl = ~2gzH(^Li + aL3) /A^ (15) 

cr = -2gz'aR/A^ (16) 

Some consequences of these operators for Higgs production in e^e~ collisions have 
been investigated in the past. Most recently, the effect of novel ZZH vertex operators 
and iiZH contact terms on the total cross sections for Higgs production has been studied 
in Ref. 0. The impact of vertex operators on angular distributions has been analyzed 
in Refs. and 0. We expand on these analyses by studying the angular distributions 
for the more general case where both novel vertex operators and contact interactions 
are present. The analysis of angular distributions in the Higgs-strahlung process ([^) 
allows us to discriminate between various novel interactions. In fact, the entire set of 
parameters ao, ai, bi and cl,cr can be determined by measuring the polar and azimuthal 
angular distributions as a function of the beam energy if the electron/positron beams 
are unpolarized. As expected, the energy dependence of the polar angular distribution is 
sufficient to provide a complete set of measurements if longitudinally polarized electron 
beams are available^. 

Total cross section and polar angular distribution. Denoting the polar angle 
between the Z boson and the e~^e~ beam axis by 6, the differential cross section for the 
process e"'"e^ r ZH may be written as 



"'^Since we can restrict ourselves to helicity-conserving couplings, as argued before, additional positron 
polarization need not be required. 
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Figure 2: Polar and azimuthal angles in the Higgs-strahlung process. [The polar angle 6^ 
is defined in the Z rest frame.] 



and the integrated cross section 



GIMI ^ A (1 + + 12 (1 + p^^^) Mils 



967rs ' " (\-Mllsf ' ' 

The Z charges of the electron are defined as usual by = — 1 and f e = — 1 + 4s^. s is the 
cm. energy squared, and A the two-particle phase space coefficient A = [1 — [jnu + rnzf I ^ 
X [1 — (rriu — mzY / s\. The coefficients a{s)^'^ and P{s)^'^ can easily be determined for 
the interactions in Eqs.(^ and (P^: 

a{sf'^ = 2ao + {s-Ml)^^c,,n (19) 

Vp ± ttp 



P(sf'R = a(s)^'^ + 27v^Mz 



AcwSw 
ai H . 

Vp ± a. 




(20) 



where the boost of the Z boson is given by 7 = (s + M| — Mfj) /2Mz^/s. 

The modification of the cross section by the new interaction terms has a simple struc- 
ture. The coefficient oq just renormalizes the SM cross section. By contrast, the contact 
interactions grow with s. [The ratio s/ h? is assumed to be small enough for the restric- 
tion to dimension-6 operators to be meaningful.] The operators O^pw, O^b affect the 
coefficient in the cross section which is independent of 6. They damp the fall-off of this 
term, changing the to a 1/ s behavior; however, these contributions remain sublead- 
ing since they are associated with transversely polarized Z bosons which are suppressed 
at high energies compared with the longitudinal components. To illustrate the size of the 
modifications a{s)^'^ and l3{s)^'^, we have depicted these functions in Fig.§(a) for the 
special choice = 1. 
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Azimuthal distributions. The azimuthal angle 0* of the fermion / is defined as the 
angle between the [e~ , Z] production plane and the [Z,f] decay plane (Fig.|^). It corre- 
sponds to the azimuthal angle of / in the Z rest frame with respect to the [e^, Z] plane. 
On general grounds, the 0* distribution must be a linear function of cos0*, cos20^,, and 
sin0*, sin 20*, measuring the helicity components of the decaying spin-1 Z state. The 
coefficients of the sine terms vanish for CP invariant theories. The cos0^, and cos 20* 
terms correspond to P-odd and P-even combinations of the fermion currents. The general 
azimuthal distributions are quite involved We therefore restrict ourselves to the 

simplified case in which all polar angles are integrated out, i.e., the polar angle 6 of the 
Z boson in the laboratory frame and the polar angle of / in the Z rest frame. In this 
way we find for the azimuthal 0* distribution: 

^ - 1 T ^ 4^ ^ (i + fr) cos 0. + ..rir^ (1 + cos 20* (21) 



32 t;2 + a2 72 + 2 V J ^* ' 2(72 + 2) 



with 



/^(7^-l)(7^-2) 



7(72 + 2) 

ai + 



4:SWCW 

Ve ± fle 


(- 


M| 

s 


Aswcw L 

Ve±ae \ 


'.") 


hi 



hi 



(22) 
(23) 



The cross section fiattens with increasing cm. energy in the Standard Model, i.e. the 
coefficients of cos0* and cos 20* decrease asymptotically proportional to and 1/s, 

respectively. The anomalous contributions modify this behavior: The cos 0* term receives 
contributions which increase proportional to ^/s with respect to the total cross section, 
while the cos 20* term receive contributions from the transversal couplings that approaches 
a constant value asymptotically. The size of the new terms in fi^^ is shown in Fig.|^(b) 
as a function of the energy. [The special choice = 1 we have adopted for illustration, 
implies f^^ = f^^.] 

High-energy limit. It is instructive to study the high-energy behavior of the coeffi- 
cients in the limit M| ^ s ^ A^. In this case we obtain the simplified relations: 

«(S)^'^ ~ TS-8swCwCL,R + 0{Ve) (24) 

/3(s)^'^ ~ a{s)^'^ + siaiT4:SwCwhi) + 0{ve) (25) 

and 

/i(s)^'^ - ^{a,T^swCwhi) + 0{ve) (26) 

/2(s)^'^ - s{a^T^swCwhl) + 0{v,) (27) 

Terms which are proportional to fg = — 1 + 4s^ are suppressed by an order of magnitude. 
If longitudinally polarized electrons are available, the asymptotic value of the coefficients 
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CLi,bi,CL and cr can be determined by measuring the polar angular distribution with- 
out varying the beam energy. The analysis of the azimuthal 0* distribution provides 
two additional independent measurements of the coefficients Oi and b^. On the other 
hand, the set of measurements remains incomplete for fixed energy if only unpolarized 
electron/positron beams are used at high energies; in this case the coefficients cannot 
be disentangled completely without varying the beam energy within the preasymptotic 
region. 
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Figure 3: Coefficients of the angular distributions as a function of the beam energy. 
Parameters are described in the text; in particular, = 1 has ben chosen in the effective 
Lagrangian Eq.^. [The L,R coefficients of the azimuthal distribution coincide for the 
special choice = 1.] 
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T. BiNOTH AND J. J. VAN DER BiJ 

Albert-Ludwigs-Universitat Freiburg, Fakultat fiir Physik, 
Hermann-Herder-Strasse 3, 79104 Freiburg, Germany 

Abstract 

We investigate the influence of massless scalar singlets on Higgs signals at the 
NLC. An exclusion bound is presented which restricts large regions of the parameter 
space but on the other hand implies that for strong interactions between the Higgs 
boson and the singlet fields of the hidden sector, detection of such a non standard 
Higgs signal can become impossible. 

1. Introduction 

Understanding of the electroweak symmetry breaking mechanism is one of the main tasks 
in particle physics. The determination of its nature would be a break-through in our 
knowledge about matter. So it is important to think about alternatives to the Standard 
Model Higgs sector. Various such extensions are available. Maybe the best motivated one 
is the supersymmetrized Standard Model with its important phenomenological implication 
of a light Higgs boson and which allows a consistent frame for grand unified theories. 
Another well understood extension - though in its minimal version disfavoured by the 
precision experiments at LEP - are technicolor theories. Though these theories avoid 
fundamental scalars, a rich bosonic spectrum of techniquark condensates may exist. Thus 
in both theories, as long as they do not occur in their minimal form, light bosonic matter 
could be present modifying the standard Higgs signals we are looking for at present and 
future colliders. If such bosons appear as singlets under the Standard Model gauge group, 
they do not feel the color or electroweak forces, but they can couple to the Higgs particle. 
As a consequence radiative corrections to weak processes are not sensitive to the presence 
of singlets in the theory, because no Feynman graphs containing singlets appear at the 
one-loop level. Since effects at the two-loop level are below the experimental precision, 
the presence of a singlet sector is not ruled out by any of the LEPl precision data. The 
only connection to such a hidden sector is a possible Higgs-singlet coupling, leading to a 
nonstandard invisible Higgs decay. The invisible decay of the Higgs boson with a narrow 
width leads to relatively sharp missing energy signals, well known from discussions on 
Majoron models [0]. However a strongly coupled hidden sector could lead to fast Higgs 
decay and thereby to wide resonances. This would disturb the signal to background ratio 
if necessary cuts are imposed. 

To check the influence of a hidden sector we will study the coupling of a Higgs boson 
to an 0(N) symmetric set of scalars, which is one of the simplest possibilities, introducing 
only a few extra parameters in the theory. The effect of the extra scalars is practically the 
presence of a possibly large invisible decay width of the Higgs particle. When the coupling 
is large enough the Higgs resonance can become wide even for a light Higgs boson. It was 
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shown earlier that there will be a range of parameters, where such a Higgs boson can be 
seen neither at LEP nor at the LHC 

In the next section we will introduce the model together with its theoretical constraints 
and in the last section we will discuss exclusion limits at the NLC 

2. The model 

The scalar sector of the model consists of the usual Higgs sector coupled to a real N- 
component vector of scalar fields, denoted by phions in the following. The Lagrangian 
density is given by, 

C = -d^<P+d^'<P-X{<i)+(P-v'^/2f-l/2d^(pdi'0-l/2m^(p^-K/(SN) {(p^f-u}/{2y/N) / 0+0 

where is the standard Higgs doublet. Couplings to fermions and vector bosons are the 
same as in the Standard Model. The ordinary Higgs field acquires the vacuum expectation 
value v/\/2. For positive oj the (^-field acquires no vacuum expectation value. After 
spontaneous symmetry breaking one is left with the ordinary Higgs boson, coupled to the 
phions into which it decays. Also the phions receive an induced mass from the spontaneous 
symmetry breaking which is suppressed by a factor 1 / If the factor N is taken to be 
large, the model can be analysed with 1/A^-expansion techniques. By taking this limit 
the phion mass remains small, but as there are many phions, the decay width of the Higgs 
boson can become large. Therefore the main effect of the presence of the phions is to give 
a large invisible decay rate to the Higgs boson. The invisible decay width is given by 

" ~ ?,2'kMh ~ 327r2aemMH 
The Higgs width is compared with the width in the Standard Model for various choices of 
the coupling uj in Fig. 0. The model is different from Majoron models 0, since the width 
is not necessarily small. The model is similar to the technicolor-like model of Ref . . 

Consistency of the model requires two conditions. One condition is the absence of a 
Landau pole below a certain scale A. The other follows from the stability of the vacuum 
up to a certain scale. An example of such limits is given in Fig. ^ where k = was taken 
at the scale 2m z-, which allows for the widest parameter range. The regions of validity 
up to a given scale A are sandwiched between the lower-left and the upper-right contour 
lines in the figure. The first stem from instability of the vacuum, the second from the 
presence of a Landau pole at that scale. 

To search for the Higgs boson there are basically two channels, one is the standard 
decay, which is reduced in branching ratio due to the decay into phions. The other is the 
invisible decay, which rapidly becomes dominant, eventually making the Higgs resonance 
wide (see Fig. |l|). In order to give the bounds we neglect the coupling n as this is a small 
effect. We also neglect the phion mass. For other values of the phion mass the bounds 
can be found by rescaling the decay widths with the appropriate phase space factor. Now 
we confront this two dimensional parameter space with the experimental potential of the 
NLC. 
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Figure 1: Higgs width in comparison with the Standard Model. 



3. NLC bounds 

At the NLC the upper hmits on the couphngs in the present model come essentially from 
the invisible decay, as the branching ratio into visible particles drops with increasing tp- 
Higgs coupling (cj), whereas for small uo one has to consider visible Higgs decays, too. 
Since the main source for Higgs production, the VTiy-fusion process, can not be used to 
look for invisible Higgs decay, one is in principle left with the Higgsstrahlung und ZZ- 
fusion reaction. For energies up to 500 GeV the Higgsstrahlungs cross section is dominant 
and is of comparable size to the ZZ-fusion process even if one is folding in the branching 
ratio B{Z e~^e~ , ^^^~). The possibility to tag an on-shell Z boson via a leptonic 
system which is extremely useful for the discrimination of possible backgrounds makes 
Higgsstrahlung to be the preferred production mechanism. Thus we only have considered 
reactions containing an on shell Z boson with its decay into e'^e~ or One should 

be aware that a few events from the huge WW background may survive 0], but that 
the Zui/ background is dominant after imposing the cuts defined below. Then the signal 
cross section is the well known Higgsstrahlungs cross section modified by the non standard 
Higgs width due to phion decay. With the invariant mass of the invisible phion system. 
Si, it has the form: 



We calculated the Zuu background with the standard set of graphs for Z production [ZZ- 
production, IViy-fusion and Z initial, final state radiation) by a Monte Carlo program 
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Figure 2: Theoretical limits on the parameters of the model in the uj vs. Mh plane. The 
contour lines correspond to the cutoff scales A = 10^^, 10^, 10^ and 10^ GeV. 

(see Ref. [^). To reduce the background we used the fact that the angular distribution of 
the Z-boson for the signal peaks for small values of | cos 9z \ in contrast to the background. 
Thus we imposed the cut | cos 6*^1 < 0.7. Because we assume the reconstruction of the 
on-shell Z-boson we use the kinematical relation Ez = (s + M| — SI)/{2^/s) between the 
Z energy and the invariant mass of the invisible system to define a second cut. Since the 
differential cross section da/dsj contains the Higgs resonance at s/ = Mfj, we impose the 
following condition on the Z energy: 

s + MI-{Mh + /^h? ^ s + Ml-iMH-AnY 

For the choice of Ah a. comment is in order. As long as the Higgs width is small, one is 
allowed to use small Ah, which reduces the background considerably keeping most of the 
signal events. But in the case of large (^-Higgs coupling, uj, one looses valuable events. 
To compromise between both effects we took Ah = 30 GeV. 

For the exclusion limits we assumed an integrated luminosity of 20/6^^. To define 
the 95% confidence level we used Poisson statistics similar to the description of Ref. [0]. 
The result is given in Fig. ^ One notices the somewhat reduced sensitivity for Mh — Mz 
due to a resonating Z boson in the ZZ background. For larger values of Mh the limit 
stems from the other Zvv backgrounds with W bosons in the t-channel and kinematical 
constrains. For large uj the signal ceases to dominate over the background because the 
Higgs peak is smeared out to an almost flat distribution. 
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Figure 3: Exclusion limits at the NLC due to Higgs searches. The dashed line corresponds 
to the invisible, the full line to all Higgs decay modes. 

We conclude from this analysis that the NLC can put further limits on the parameter 
space of our invisible Higgs model. Note that within the kinematic range very strong limits 
on uj can be set. Again there is a range where the Higgs boson will not be discovered, 
even if it does exist in this mass range. This has already been shown for the Higgs search 
at LEP and also holds true for the heavy Higgs search at LHC We see that a sufficiently 
wide nonstandard Higgs resonance would make it very difficult to test the mechanism of 
electroweak symmetry breaking at future colliders. 
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Abstract 

The production mechanisms and decay modes of the heavy neutral and charged 
Higgs bosons in the Minimal Supersymmetric Standard Model are investigated 
at future e^e~ colliders in the TeV energy regime. We generate supersymmet- 
ric particle spectra by requiring the MSSM Higgs potential to produce correct 
radiative electroweak symmetry breaking, and we assume a common scalar 
mass mo, gaugino mass mi/2 and trilinear coupling A, as well as gauge and 
Yukawa coupling unification at the Grand Unification scale. Particular em- 
phasis is put on the low tg/3 solution in this scenario where decays of the Higgs 
bosons to Standard Model particles compete with decays to supersymmetric 
charginos/neutralinos as well as sfermions. In the high tg/3 case, the super- 
symmetric spectrum is either too heavy or the supersymmetric decay modes 
are suppressed, since the Higgs bosons decay almost exclusively into h and r 
pairs. The main production mechanisms for the heavy Higgs particles are the 
associated AH production and H^H~ pair production with cross sections of 
the order of a few fb. 



1. Introduction 

Supersymmetric theories are generally considered to be the most natural extensions 
of the Standard Model (SM). This proposition is based on several points. In these theo- 
ries, fundamental scalar Higgs bosons with low masses can be retained in the context 
of high unification scales. Moreover, the prediction of the renormalized electroweak 
mixing angle sin^ 9w = 0.2336 ± 0.0017, based on the spectrum of the Minimal Super- 
symmetric Standard Model (MSSM) |T0[, is in striking agreement with the electroweak 
precision data which yield sin^ 9w = 0.2314(3) |TI|. An additional attractive feature is 
provided by the opportunity to generate the electroweak symmetry breaking radiatively 
[^. If the top quark mass is in the range between ~ 150 and ~ 200 GeV, the universal 
squared Higgs mass parameter at the unification scale decreases with decreasing energy 
and becomes negative at the electroweak scale, thereby breaking the SU(2)l x U(1)y 
gauge symmetry while leaving the U(l) electromagnetic and SU(3) color gauge symme- 
tries intact 0. The analysis of the electroweak data prefers a light Higgs mass |1TT|,|9[| as 
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predicted in supersymmetric theories; however since the radiative corrections depend only 



logarithmically on the Higgs mass |T^, the dependence is weak and no firm conclusions 
can yet be drawn. 

The more than doubling the spectrum of states in the MSSM gives rise to a rather 
large proliferation of parameters. This number of parameters is however reduced drasti- 
cally by embedding the low-energy supersymmetric theory into a grand unified (GUT) 
framework. This can be achieved in supergravity models in which the effective low- 
energy supersymmetric theory [including the interactions which break supersymmetry] is 
described by the following parameters: the common scalar mass mo, the common gaugino 
mass mi/2, the trilinear coupling A, the bilinear coupling B, and the Higgs-higgsino mass 
parameter /i. In addition, two parameters are needed to describe the Higgs sector: one 
Higgs mass parameter [in general the mass of the pseudoscalar Higgs boson. Ma] and the 
ratio of the vacuum expectation values, tg/3 = f2/fi, of the two Higgs doublet fields which 
break the electroweak symmetry. 

The number of parameters can be further reduced by introducing additional con- 
straints which are based on physically rather natural assumptions: 



(z) Unification of the b and r Yukawa couplings at the GUT scale [11| leads to a 



correlation between the top quark mass and tg/3 [^|13],0]. Adopting the central value 



of the top mass as measured at the Tevatron [|1^ , tg(3 is restricted to two narrow ranges 



around tg/3 ~ 1.7 and 50, with the low tg(3 solution theoretically somewhat favored |jT4 

[a) If the electroweak symmetry is broken radiatively, then the bilinear coupling B 
and the Higgs-higgsino mass parameter n are determined up to the sign of /i. [The 
sign of /i might be determined by future precision measurements of the radiative b decay 
amplitude.] 

{iii) It turns out a posteriori that the physical observables are nearly independent of 
the GUT scale value of the trilinear coupling Aq, for \Ag\ < 500 GeV. 

Mass spectra and couplings of all supersymmetric particles and Higgs bosons are 
determined after these steps by just two mass parameters along with the sign of fi; we 
shall choose to express our results in terms of the pseudoscalar Higgs boson A mass Ma 
and the common GUT gaugino mass mi/2- 

In this paper we focus on heavy Higgs particles A, H and i?^ with masses of a few 
hundred GeV, and therefore close to the decoupling limit ||16[. The pattern of Higgs 
masses is quite regular in this limit. While the upper limit on the mass of the lightest 
CP-even Higgs boson /i is a function of tg/3 [jl^ , 



Mh < 100 to 150 GeV [for low to high tg^] (1.1) 
the heavy Higgs bosons are nearly mass degenerate [c.f. Fig.l] 

Ma ^ Mh ^ Mh± (1.2) 
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Moreover, the properties of the hghtcst CP-even Higgs boson h become SM-hke in this 
hmit. The production of the heavy Higgs bosons becomes particularly simple in e'^e~ 
collisions; the heavy Higgs bosons can only be pair-produced, 



e+e" ^ AH (1.3) 
e+e" ^ H+H- (1.4) 

Close to this decoupling hmit, the cross section for H Higgs-strahlung e+e~ — > ZH is very 
small and the cross section for the WW fusion mechanism e'^e~ — > fe^^eH is appreciable 



only for small values of tg/?, tg/? ~ 1, and relatively small H masses, Mh ^ 350 CcV. The 
cross section for ZZ fusion of the H is suppressed by an order of magnitude compared to 
WW fusion. The pseudoscalar A particle does not couple to W/Z boson pairs at the tree 
level. 

The decay pattern for heavy Higgs bosons is rather complicated in general. For large 
tg/3 the SM fermion decays prevail. For small tg/3 this is true above the tt threshold of 
Mh,a <i 350 GeV for the neutral Higgs bosons and above the tb threshold of Mh± ^180 
GeV for the charged Higgs particles. Below these mass values many decay channels 
compete with each other: decays to SM fcrmions // [and for H to gauge bosons VV], 
Higgs cascade decays, chargino/neutralino XiXj decays and decays to supersymmetric 
sfermions // 

H ^ ff,VV,hh, XiXj , /? (1-5) 
A ^ ff,hZ, XiXj , ff (1.6) 
^ //', hW^ , XiXj, ff' (1.7) 

In this paper, wc analyze in detail the decay modes of the heavy Higgs particles and 
their production at e+e^ linear colliders. The analysis will focus on heavy particles for 
which machines in the TeV energy range are needed. The paper is organized in the 
following way. In the next section we define the physical set-up of our analysis in the 
framework of the MSSM embedded into a minimal supergravity theory. In section 3, 
we discuss the production cross sections of the heavy Higgs bosons. In the subsequent 
sections, we discuss the widths of the various decay channels and the final Higgs decay 
products. 



2. The Physical Set-Up 

The Higgs sector of the Minimal Supersymmetric Standard Model is described at tree-level 
by the following potential 



Vo = (mjj^+ii')\H,\^ + (mjj^+ii'')\H2\^-ml(eijH^'H2^ + h.c.) 

1 r n2 1 



+ 2^ i^r^2 



i|2 



(2.1) 
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The quadratic Higgs terms associated with /i and the quartic Higgs terms coming with 
the electroweak gauge couphngs g and g' are invariant under supersymmetric transforma- 
tions, m^^, ''^H2 "^3 soft-supersymmetry breaking parameters with mg = Bfi. 
eij [i,j = 1,2 and ei2 = 1] is the antisymmetric tensor in two dimensions and Hi = 
[Hi, Hi) = {Hi, Hi), H2 = {Hi Hi) = {H^,H^) are the two Higgs-doublet fields. After 
the symmetry breaking, three out of the initially eight degrees of freedom will be absorbed 
to generate the and Z masses, leaving a quintet of scalar Higgs particles: two CP- 
even Higgs bosons h and H, a CP-odd [pseudoscalar] boson A and two charged Higgs 
particles H"^. 

Retaining only the [leading] Yukawa couplings of the third generation 

Xt = — — r , Afe = and = (2.2) 

f sm p V cos p V cos p 

where tg/5 = V2/V1 [with = vf + fixed by the W mass, v = 246 GeV] is the ratio 
of the vacuum expectation values of the fields H2 and H^, the superpotential is given in 
terms of the superfields Q = {t, b) and L = (r, z/,-) byQ 



W = e,,- 



XtHiQH" + XbHiQ^W + XrH[UT^ - fiHlHi (2.3) 



Supersymmetry is broken by introducing the soft-supersjTiimetry breaking bino B, wino 
[a =1-3] and gluino g"' [a =1-8] mass terms, 

^MiBB + ^M^WW^ + ^M,Tg\ (2-4) 

soft-supersymmetry breaking trilinear couplings, 

eij [XtAtHiQH^ + XbAbHlQ^b' + XrArH\Vr - fiBHlHi] (2.5) 

and soft-supersymmetry breaking squark and slepton mass terms 

[ilh + blbL] + Mlji*^tR + MlVj^bR + Ml [fl^L + i^rli^ri] + Mlr^fn + ■■■ (2.6) 

where the ellipses stand for the soft mass terms corresponding to the first and second 
generation sfermions. 

The minimal SUSY-GUT model emerges by requiring at the GUT scale Mq'- 

{i) the unification of the U(l), SU(2) and SU(3) couphng constants = gf/4:n [i = 1-3], 

a3(MG) = a2{MG) = ai{MG) = aa (2.7) 



^Note that our convention for the sign of /i is consistent with Ref.[|l^, which is opposite to the one 
adopted in Ref.[|9|. 
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(ii) a common gaugino mass; the Mi with i =1-3 at the electroweak scale are then related 
through renormalization group equations (RGEs) to the gauge couplings, 

M. = '^^m,,, _ AUAh) = ^M,(M.) = '^^MdAh) (2.8) 

OiG a2{Mz) ai{Mz) 

{Hi) a universal trilinear coupling A 

Ac = At{Mo) = A(Mg) = A.(Mg) (2.9) 

[Hi) a universal scalar mass mo 

mo = Mq = Mu = Md = Ml = Me 

= muAMG) = mH,{MG) (2.10) 

Besides the three parameters mi/2,AG and mo the supersymmetric sector is described 
at the GUT scale by the bilinear coupling Bq and the Higgs-higgsino mass parameter 
liG- The theoretically attractive assumption that the electroweak symmetry is broken 
radiatively constrains the latter two parameters. Indeed, radiative electroweak symmetry 
breaking results in two minimization conditions [see Ref.|jT9| for details] of the Higgs 
potential; at the low-energy scale in the tree approximation, they are given by 

-Ml = ^ --/i^ 2.11 

2 ^ tan2/?-l ^ ^ ' 

Bfi = ^{mjj^+mjj^ + 2fi^)sm2p (2.12) 

For given values of the GUT parameters mi/2, mo, Ac as well as tg/3, the first minimization 
equation can be solved for /i [to within a sign]; the second equation can then be solved 
for B. Since m|^^ and m^^ are related to Ma through the RGEs, the solution for n 
and B can be approximately expressed as a function of Ma and tg/3. The power of 
supergravity models with radiative electroweak symmetry breaking becomes apparent 
when one includes the one-loop contributions to the Higgs potential. It is through these 
one-loop terms that most of the supersymmetric particle masses are determined; the 
minimization conditions [which are also solved for fi to within a sign and B] fix the 
masses in order that the electroweak symmetry is broken correctly, i.e. with the correct 
value of Mz- [U(1)em and SU(3) remain unbroken of course]. The one-loop contributions 
and the minimization equations are given in Ref. to which we refer for details. 



A heavy top quark is required to break the electroweak symmetry radiatively, since 
it is the large top Yukawa coupling which will drive one of the Higgs mass parameters 
squared to a negative value. As emphasized before, the additional condition of unification 
of the b-T Yukawa couplings gives rise to stringent constraints on tg(3. The attractive 
idea of explaining the large top Yukawa coupling as a result of a fixed point solution of 
the RGEs leads to a relationship between Mt and the angle jS, Mt ~ (200 GeV) sin /3 for 
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tg/5 < 10, giving a further constraint on the model. 

To hmit the parameter space further, one could require that the SUGRA model is 
not fine-tuned and the SUSY breaking scale should not be too high, a constraint which 
can be particularly restrictive in the small tg(3 region. However, the degree of fine-tuning 
which can be considered acceptable is largely a matter of taste, so we disregard this issue 
in our analysis. 

We now detail the calculations of the supersymmetric particle spectrum more precisely. 
We incorporate boundary conditions at both electroweak and GUT scales, adopting the 
ambidextrous approach of Ref . |jT9[ . We specify the values of the gauge and Yukawa cou- 



plings at the electroweak scale, in particular M^, tg(3 and a^. The gauge and Yukawa 
couplings are then evolved to the GUT scale Mq [defined to be the scale jl for which 



Q!i(/i) = a2{jl)] using the two-loop RGEs [|1^]. At Mq we specify the soft supersym- 
metry breaking parameters mi/2, rriQ and Aq. We then evolve parameters down to the 
electroweak scale where we apply the one-loop minimization conditions derived from the 
one-loop effective Higgs potential and solve for /i to within a sign and B [we then can 
integrate the RGEs back to Mq and obtain /ic and Bq]- By this procedure, the su- 
persymmetric spectrum is completely specified; that is, we generate a unique spectrum 
corresponding to particular values of mi/2, mo, Aq, tg/5 and the sign of fi. It turns out 
that the spectrum is nearly independent of Aq, for \Ag\ < 500 GeV. In most of our cal- 
culations, we substitute a particular value of Ma for mo in order to introduce a mass 
parameter which can be measured directly. 

We discuss the SUSY spectrum and its phenomenological implications for two repre- 
sentative points in the M^-tan/? planeQ. We choose Mf"^"^ = 175 GeV, consistent with the 
most recent Tevatron analyses ||T5| throughout our calculations, and values of tg/3 = 1.75 
and 50, which are required (within uncertainties) by 6-r unification at Mq. In particular, 
we emphasize the low tg(3 solutions; they are theoretically favored from considerations 
such as b —* s'y [^ and cosmological constraints [^. The low tg(3 solutions generate 



much lighter SUSY spectra, more likely to be seen at future e~^e~ colliders. In both 
the low and high tg/3 regions we take^ ^^(M^) = 0.118 [|^ and Aq = 0, though the 
qualitative behavior in each region should not depend greatly on these parameters, 
(a) Low tan f3 

As a typical example of the low tan /3 region we consider the point Mf = 175 GeV and 



tan/5 = 1.75 for which Xt{MG) lies in its "fixed-point" region [p!2| , p!^ . If Ma is fixed, the 
scalar mass parameter mo can be calculated as a function of the common gaugino mass 
parameter mi/2 so that all Higgs and supersymmetric particle masses can in principle be 
parameterized by mi/2- The correlation between mo and mi/2 is shown in Fig. 2 for three 

^Our numerical analysis is consistent with the numbers obtained in Ref. [pO| , once their value of A-r in 
Tab. 2 is corrected. We thank W. de Boer for a discussion on this point. 

•^This corresponds to the sin d\Y value quoted and compared with the high-precision electroweak 
analyses in the Introduction. 
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values of Ma = 300, 600 and 900 GeV in the low igjS region. 

Some of the parameter space is already eliminated by experimental bounds on the 
light Higgs mass, the chargino/neutralino masses, the light stop mass, the slepton masses 



and the squark/gluino masses from LEPl/1.5 and the Tevatron [|2^. The lower limits are 
indicated by the non-solid lines in Fig. 2. Low values of mi/2 ^ 60 GeV are excluded by 
the lower bound on the gaugino masses. For > 0, the bound from the negative search 
of charginos at LEP1.5 almost rules out completely the scenario with Ma ^ 300 GeV. If 
the h boson is not discovered at LEP2, i.e. if M^ > 95 GeV, the whole /i < scenario [for 
^1/2, '"T'O < 500 GeV] can be excluded, while for /x > only the mi/2 > 200 GeV range 
[which implies very large values of Ma] would survive. The requirement that the lightest 
neutralino is the LSP, and therefore its mass is larger than the lightest f mass, excludes 
a small edge of the parameter space [dotted line] at small mo with mi/2 > 200 GeV in 
the /i < case. 

The masses of the Higgs bosons are shown in Fig. 3a as a function of mi/2 for tgf3 = 
1.75, both signs of fi and for two representative values of mo = 100 and 500 GeV. The 
lightest Higgs boson has a rather small tree-level mass and M^ comes mainly from ra- 
diative corrections; the maximal values [for mi/2 ~ 400 GeV] are M^'^^ ~ 90 GeV for 
/i < and ~ 100 for /i > 0. Because the pseudoscalar mass is approximately given 
by ~ Bfi/ sin 2/3 ~ Bfi [at the tree-level] and since Bfi turns out to be large 
in this scenario, the pseudoscalar A is rather heavy especially for large values of mo, 
and thus is almost mass degenerate with the heavy CP-even and charged Higgs bosons. 
Ma ~ ~ Mh±. Note that Ma is below the ti threshold. Ma < 350 GeV, only if mo 
and mi/2 are both of (9(100) GeV. 

The chargino/neutralino and sfermion masses are shown Fig.3b-d as a function of 
mi/2 for the two values Ma = 300 and 600 GeV and for both signs of fi. In the case of 
charginos and neutralinos, the masses are related through RGEs by the same ratios that 
describe the gauge couplings at the electroweak scale. The LSP is almost bino-like [with 
a mass m^o ~ Mi] while the next-to-lightest neutralino and the lightest chargino are 
wino-like [with masses m^o ~ m^+ ~ M2 ~ 2m^o]. The heavier neutralinos and chargino 
are primarily higgsinos with masses m^o ~ m^o ~ m^+ ~ Note that the masses 
approximately scale as Ma and that the decay of the heavy scalar and pseudoscalar Higgs 
bosons into pairs of the heaviest charginos and neutralinos is kinematically not allowed. 

The left- and right-handed charged sleptons and sneutrinos are almost mass degen- 
erate, the mass differences not exceeding (9(10) GeV; the mixing in the r sector is rather 
small for small tgjS, allowing one to treat all three generations of sleptons on the same 
footing. In the case of squarks, only the first two generations are degenerate, with left- 
and right-handed squarks having approximately the same mass. The mixing in the stop 
as well as in the sbottom sector leads to a rather substantial splitting between the two 
stop or sbottom mass eigenstates. Only for small values of Ma and for /i < is bi the 
lightest squark; otherwise ii is the lightest squark state. Note that the squark masses 
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increase with mi/2 and that they scale as Ma i.e. as The slepton masses decrease 
with increasing mi/2: this is due to the fact that when mi/2 increases and Ma is held 
constant, uiq decreases [see Fig. 2], and the dependence of the slepton masses on mo is 
stronger [for fixed mo, the slepton masses would increase with increasing mi/2]. 
(b) High tan (5 

In this region we take tg/3 = 50 as a representative example, a value consistent with 
the unification of the t, h and r Yukawa couplings. The set of possible solutions in the 
parameter space [mi/2, "'^o] for = 300 and 600 GeV is shown in Fig. 4. At tg/3 = 50 and 
^poie _ j^-j^g GeV, we find solutions only for /i < 0; this is a result of the large one-loop 
contribution to M^, the sign of which depends on [^. The boundary contours given 
in the figure correspond to tachyonic solutions in the parameter space: m^^ < 0, M\ < 
or < at the tree-level. The latter constraint is important for algorithmic reasons: 
M| at the tree-level enters into the minimization equations in the form of a logarithm 



1^. Also the requirement of the lightest neutralino to be the LSP excludes a small edge 



of the parameter space at small values of mo; this explains why the curves for M4 = 300 
and 600 GeV do not terminate for low mo values. 



Particle 


Mass (GeV) 


Mass (GeV) 


Mass (GeV) 


Mass (GeV) 


Ma 


300 


300 


600 


600 


(mi/2, mo) 


(364,250) 


(352,800) 


(603,300) 


(590,800) 


9 


940 


910 


1557 


1524 




662,817 


753,896 


1115,1285 


1156,1325 




689,787 


804,894 


1159,1260 


1220,1312 


Ul,U2 


881,909 


1144,1164 


1431,1479 


1586,1628 




878,912 


1142,1167 


1425,1481 


1582,1630 


n, T2; 


165,365; 325 


567,740; 729 


255,517; 485 


586,812; 799 


ei, 62; l?e 


290,360; 351 


813,838; 835 


381,519; 513 


833,901; 898 


xt 


268,498 


261,536 


452,764 


443,779 




144,268,485,496 


139,261,526,534 


239,452,754,763 


234,443,771,778 


MA,MH±,MH,Mh 


300,315,300,124 


300,315,300,124 


600,608,600,131 


600,608,600,131 



Tab.l: Particle spectra for Mf°'^ = 175 GeV, tan/3 = 50 for selected M^,mi/2 aiid tuq values. 



The sparticle spectra for Ma = 300 and 600 GeV and two sets of mi/2 and (extreme) 
mo values are shown in Table 1. In all these cases, the particle spectrum is very heavy; 
hence most of the SUSY decay channels of the Higgs particles are shut for large tg/3. The 
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only allowed decay channels are H,A-^ fifi, XiXi and H"^ fiu [for large Ma values]. 
However, the branching ratios of these decay channels are suppressed by large hh and 
th widths of the Higgs particles for large tg/3: while the supersymmetric decay widths 
are of the order (9(0.1 GeV), the decays involving h quarks have widths (9(10 GeV) and 
dominate by 2 orders of magnitude. 

(c) Additional Constraints 

There are additional experimental constraints on the parameter space for both high and 
low tg/3; the most important of these are the h — >■ S7, Z bb, and dark matter [relic LSP 
abundance] constraints. These constraints are much more restrictive in the high tg/3 case. 



Recent studies [BT| have indicated that the combination of 6 ^ 57, dark matter and 



nih constraints disfavor the high tg/9 solution for which the t, b and r Yukawa couplings are 
equal, in particular the minimal SUSY-SO(IO) model with universal soft-supersymmetry 
breaking terms at Mq- This model can, however, be saved if the soft terms are not univer- 
sal [implying a higgsino-like lightest neutralino], and there exist theoretical motivations 
for non-universal soft terms at Mq [^. The presently favored Z ^ bb decay width would 
favor a very low A mass for large tg/3. 

For low tg/3, these additional constraints do not endanger the model, yet they can 
significantly reduce the available parameter space. In particular the Z — > 66 constraint 
favors a light chargino and light stop for small to moderate values of tg/3 [^^ so that 
they could be detected at LEP2 The dark matter constraint essentially places an 

upper limit on itlq and mi/2 ||2^. The b s'j constraint [^, on the other hand, is 



plagued with large theoretical uncertainties mainly stemming from the unknown next-to- 
leading QCD corrections and uncertainties in the measurement of as{Mz)- However, it is 
consistent with the low tg(3 solution and may in the future be useful in determining the 



sign of /i [pi . 



3. Production Mechanisms 

The main production mechanisms of neutral Higgs bosons at e~^e~ colliders are the Higgs- 
strahlung process and pair production, 

(a) Higgs-strahlung e^e^ {Z)^Z + h/H 
(6) pair production e^e^ {Z) ^ A + h/H 

as well as the WW and ZZ fusion processes, 

(c) fusion processes e'^e~ uv [WW) ^ uu -\-h/H 

e+e" ^ e+e-{ZZ) e+e" + h/H 

[The CP-odd Higgs boson A cannot be produced in the Higgs-strahlung and fusion pro- 
cesses to leading order since it does not couple to VV pairs.] The charged Higgs particle 
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can be pair produced through virtual photon and Z boson exchange, 



(d) charged Higgs e^e (7, Z*) H^H 

[For masses smaller than ~ 170 GeV, the charged Higgs boson is also accessible in top 
decays, e~^e~ tt with t H^h.] 

The production cross sections^ for the neutral Higgs bosons are suppressed by mixing 
angle factors compared to the SM Higgs production, 

cr(e+e- ^ Zh) , aiVV h) , cr(e+e- ^ AH) ~ sm^{p - a) (3.1) 
a(e+e- ^ ZH) , a{yV H) , a^e+e' Ah) ~ cos^(/3 - a) (3.2) 

while the cross section for the charged Higgs particle does not depend on any parameter 
other than Mh±. 

In the decoupling limit. Ma ^ M^, the HVV couphng vanishes, while the hVV 
coupling approaches the SM Higgs value 

gjjvv = cos{(3 -a) ^ Ml sin 4/3/2M^ ^ (3.3) 
ghvv = sm{(3-a)^l-0{Ml/M^) ^1 (3.4) 

Hence, the only relevant mechanisms for the production of the heavy Higgs bosons in this 
limit will be the associated pair production (b) and the pair production of the charged 
Higgs particles (d). The cross sections, in the decouphng limit and for ^ Mz, are 
given by [we use Mh ~ Ma] 

a{e^e-^AH) = ^{v^ + al)Pl (3.5) 



a(e+e-^H+H-) = ^^^^^^^^w 



^ ^ o 2 2 ^ 9e;R^4 4 



(3.6) 



where /3j = (1 — 4M|/s)^/^ is the velocity of Higgs bosons, the Z couplings to electrons 
are given by Og = — l,fe = — 1 + 4sin^6'i4/, and to the charged Higgs boson by vu = 
-2 + 4sin2^VK- The cross sections for hA and HZ production vanish in the decoupling 
limit since they are proportional to cos^(/3 — a). 

The cross section for the fusion process, e~^e~ — >■ Uei^eH, is enhanced at high energies 
since it scales like M^ \og s/Mfj. This mechanism provides therefore a useful channel for 
H production in the mass range of a few hundred GeV below the decoupling limit and 
small values of tg/?, where cos^(/9 — a) is not prohibitively small; the cross section, though, 

^The complete analytical expressions of the cross sections can be found, e.g., in Ref.[^. Note that in 
this paper there are a few typos that we correct here: in eq.(20), the factor 92 should replaced by 96; in 
the argument of the A function of the denominator in eq.(21), the parameter M\ should be replaced by 
M"^-. finally, the minus sign in the interference term in eq.(25) should be replaced by a plus sign. 
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becomes gradually less important for increasing Mh and vanishes in the decoupling limit. 
In the high energy regime, the WW — *• H fusion cross section is well approximated by 
the expression 



/ + - - TT\ G'pMy^f 



3 



' M2\ s ( Ml' 



coi?{(3-a) (3.7) 



obtained in the effective longitudinal W approximation. Since the NC couplings are small 
compared to the CC couplings, the cross section for the ZZ fusion process is ~ 16 cos"^ 6w, 
i.e. one order of magnitude smaller than for WW fusion. 

Numerical results for the cross sections are shown in Fig.5 at high-energy e^e~ colliders 
as a function of ^/s TeV for the two values tg/5 = 1.75 and 50, and for pseudoscalar masses 
Ma = 300, 600 and 900 GeV [note that Mh ^ Mh± ^ Ma as evident from Figs. 1 and 
3a]. For a luminosity oi J C = 200 fb~^, typically a sample of about 1000 HA and H^H~ 
pairs are predicted for heavy Higgs masses of ~ 500 GeV at ^/s = 1.5 TeV. For small tg/3 
values, tg/3 < 2, a few hundred events are predicted in the WW — > H fusion process for 
H masses ~ 300 GeV. The cross sections for the hA and HZ processes are too low, less 
than ~ 0.1 fb, to be useful for Mh > 300 GeV; Fig.5b. 

Note that the cross sections for the production of the lightest Higgs boson h in the 
decoupling limit and for ^ Mz, Mh are simply given by 

aie^e-^ZZ) ^ ^(-e^ + ^D (3-8) 



G%M, 



s 



^ ' 4^2773 ^ Ml ^ ' 

The cross sections are the same as for the SM Higgs particle and are very large ~ 100 fb, 
especially for the WW fusion mechanism. 



4. Decay Modes 

4.1 Decays to standard particles 

For large tg(3 the Higgs couplings to down-type fermions dominate over all other 
couplings. As a result, the decay pattern is in general very simple. The neutral Higgs 
bosons will decay into bb and r+r^ pairs for which the branching ratios are close to ~ 90 % 
and ~ 10 %, respectively. The charged Higgs particles decay into rur pairs below and 
into tb pairs above the top-bottom threshold. 

The partial decay widths of the neutral Higgs bosons[|, $ = if and A, to fermions are 
given by 

r(t - //) = N^^djfnjir, (4.1) 

^We refrain from a discussion of the h decays which become SM-like in the decoupUng hmit. In 
addition, we discuss only the dominant two-body decay modes of the heavy Higgs bosons; for an updated 
and more detailed discussion, including also three-body decays, see Ref . p3| . 
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with p = 3(1) for the CP-even (odd) Higgs bosons; /?/ = (1 — 4mj/M|)-^/^ is the velocity 
of the fermions in the final state, Nc the color factor. For the decay widths to quark pairs, 
the QCD radiative corrections are large and must be included; for a recent update and a 
more detailed discussion, see Ref . [01 . 



The couplings of the MSSM neutral Higgs bosons [normalized to the SM Higgs coupling 



9HsMff 



r- 1/2 

\/2Gp ruf and gng 



1/2 



are given in Table 2. 





Q^uu 


9^dd 


g^vv 


h 


COS a/ sin /5 


— sin a/ cos [3 


sin(/5 — a) 


H 


sin a/ sin /5 


cos a/ cos/3 


cos(/3 — a) 


A 


l/tg/3 








Tab. 2: Higgs boson couplings in the MSSM to fermions and gauge bosons relative to the SM 
Higgs couplings. 



In the decoupling limit. Ma ^ M^, we have 



cos a ~ sin/3 — cos/3— — 2-sin4/9 — sin/9 



2Ml 



sin a ~ — cos /3 + sin /3 „ sin 4/3 — » — cos/3 



2Ml 



(4.2) 
(4.3) 



Therefore the hf f couplings reduce to the SM Higgs couplings, while the Hff couplings 
become equal to those of the pseudoscalar boson A, 



cos a/ sin (3 

— sin a/ cos /3 

— sin a I sin /3 
cos a I cos /3 

The partial width of the decay mode 



1 
1 

- l/tg/3 

- tg/3 

Mc? is given by 



(4.4) 



T{H+ ^ ud) 



4^271 Mj^± 
"(M 



) (m^tg'/3 + m^ctg2/3) - 4m2 



(4.5) 



with Vurf the CKM-type matrix element for quarks and A is the two-body phase space 
function defined by 



A,,, = (1 - Ml /Ml - Mf/Mlf - AMlMf/M^ 
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(4.6) 



For decays into quark pairs, the QCD corrections must be also included. 

Below the tt threshold, a variety of channels is open for the decays of the heavy CP- 
even Higgs bosons, the most important being the cascade decays H — > with $ = /i or 
A, with a partial width [for real light Higgs bosons] 

where = (1 — 4M|/M^)^/^ and the radiatively corrected three-boson self-couplings 
[to leading order], in units of g'^ = {y/^GpY^'^M'^, are given by 

2 

6 sm (X cos cx 

9Hhh — 2sin2Q;sin(/3 + a) — cos2q;cos(/3 + a) + 3— J- : — ^ (4.8) 



Qhaa = - cos 2/3 cos(/9 + a) + 



Mi sin (3 

e sin a cos^ (3 



M| sin f3 

In contrast to the previous couplings, the leading radiative corrections cannot be 
absorbed entirely in the redefinition of the mixing angle a, but they are renormalized by 
an explicit term depending on the parameter e given by [Ms is the common squark mass 
at the electroweak scale] 



^=^-T^log 1 + (4-9) 



V^7r2 sin^ P \ m, 



In the decoupling limit, these couplings approach the values 

3 

gnhh 2 

9HAA sin 4/3 (4.10) 

In the mass range above the WW and ZZ thresholds, where the HVV couplings are 
not strongly suppressed for small values of tg/3, the partial widths of the H particle into 
massive gauge bosons can also be substantial; they are given by 

T{H VV) = ^^^'"f^"~^V |,(l - Akv + 124)(1 - AKyf' 5'y (4.11) 

with Kv = M'^/Ml and 5'y = 2(1) for V = W{Z). 

For small values of tg/3 and below the tt and the tb thresholds, the pseudoscalar and 
charged Higgs bosons can decay into the lightest Higgs boson h and a gauge boson; 
however these decays are suppressed by cos^(/3 — a) and therefore are very rare for large 
A masses. The partial decay widths are given by 
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In the decoupling limit, the partial widths of all decays of the heavy CP-even, CP- 



odd and charged Higgs bosons involving gauge bosons vanish since cos^(/3 



a) 



0. In 



addition, the H ^ hh decay width is very small since it is inversely proportional to Mh, 
and H AA is not allowed kinematically. Therefore, the only relevant channels are the 
decays into bb/tt for the neutral and tb for the charged Higgs bosons. The total decay 
widths of the three bosons H, A and H^, into standard particles can be approximated in 
this limit by 

3Gp 



T{Hk ^ all) 



4v^7r 



(4.13) 



[We have neglected the small contribution of the decays into r leptons for large tgjS. 



4.2 Decays to charginos and neutralinos 



The decay widths of the Higgs bosons [k = (1, 2, 3, 4) correspond to (if, h, A, H^) 
into neutralino and chargino pairs are given by 

,2 \ 



r(i/fc XiXj) 



GfM^Mh.X. 



1/2 



2^71 1 + 5i 



'^Vk^i^j Pij k Fj i k 



m 



Xi 



m 



Ml 



(4.14) 



where ?7i,2,4 = +1, ^73 = —1 and 5ij = unless the final state consists of two identical 
(Majorana) neutralinos in which case da = 1; = ±1 stands for the sign of the i'th 
eigenvalue of the neutralino mass matrix [the matrix Z is defined in the convention of 
Ref. |jl8|, and the eigenvalues of the mass matrix can be either positive or negative] while 
ej = 1 for charginos; Xij^k is the usual two-body phase space function given in eq.(4.4). 

In the case of neutral Higgs boson decays, the coefficients Fijk are related to the 
elements of the matrices U, V for charginos and Z for neutralinos. 



Xi Xj 



Xi Xj 



Fijk — [ekViiUj2 — dkVi2Uji] 

Fijk = - {Zj2 — tanOwZji) {ckZi^ + dkZii) + i 



(4.15) 
(4.16) 



with the coefficients and given by 

ei/(ii = cosa/ — sina , e2/c?2 = sina/cosa , 63/^3 = — sin /5/ cos/? (4-17) 
For the charged Higgs boson, the coupling to neutralino/chargino pairs are given by 



jii 



cos P 
sin P 



ZjiVii + —j= {Zj2 



ianOwZji) Vi 



i2 



ZjsUii — — ^ {Zj2 + tanOwZji) Ui2 



(4.18) 
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The matrices U, V for charginos and Z for neutralinos can be found in Ref . |TB . 



Since in most of the parameter space discussed in Section 2, the Higgs-higgsino mass 
parameter |/i| turned out to be very large, 3> Mi, M2, Mz, it is worth discussing the 
Higgs decay widths into charginos and neutrahnos in this hmit. First, the decays of the 
neutral Higgs bosons into pairs of [identical] neutralinos and charginos if^ —>■ XiXi will 
be suppressed by powers of M|//i^. This is due to the fact that neutral Higgs bosons 
mainly couple to mixtures of higgsino and gaugino components, and in the large /i limit, 
neutralinos and charginos are either pure higgsino- or pure gaugino-like. For the same 
reason, decays — > Xi 2X1 and Xs 4X2 of the charged Higgs bosons will be suppressed. 
Furthermore, since in this case Ma is of the same order as decays into pairs of heavy 
charginos and neutralinos will be kinematically forbidden. Therefore, the channels 

H,A X? X3,4 ' xl X3,4 and xt X2 

^ Xtxl, andx^X?,2 (4-19) 

will be the dominant decay channels of the heavy Higgs particles. Up to the phase space 
suppression [i.e. for M^ sufficiently larger than the partial widths of these decay 
channels, in units of Gi^M^Mj^^ /(4v^7r), are given by ||35[ 

V{H^xIxIa) = ^tan2^^(l±sin2/5) 

m^xlxl,) = ^(1± sin 2/3) 

V{H^xfx^) = 1 (4.20) 

TiA^XiXh) = ^tan2^H^(l±sin2/?) 

m^xlxh) = ^(1± sin 2/5) 

T{A^XiX^) = 1 (4.21) 

m^^xtxh) = 1 

r{H+^xtxi) = 1 

TiH+^xtxl) = tan'Ow (4.22) 

[We have used the fact that in the decoupling limit sin 2a = — sin 2(3.] If all these channels 
are kinematically allowed, the total decay widths of the heavy Higgs bosons to chargino 
and neutralino pairs will be given by the expression 

3GpM^^^ / , 1 . 



TiHk ^ E x^XJ) = ^^J^ Mh, + - tan^ 9^) (4.23) 

which holds universally for all the three Higgs bosons H, A, H^. 
4.3 Decays to squarks and sleptons 
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Decays of the neutral and charged Higgs bosons, Hk = h, H, A, H^, to sfermion pairs can 
be written as 



T{Hk — > fifj 



NcGp yl/2 



(4.24) 



fi with 2 = 1,2 are the sfermion mass eigenstates which are related to the current eigen- 
states /l, Jr by 



/i = fLCOsOf + fRsinOf 
h = -hsmOf + fRCOs9f 



(4.25) 



The mixing angles Of are proportional to the masses of the partner fermions and they 
are important only in the case of third generation sfermions. The couplings Quf^f j' 
neutral and charged Higgs bosons to sfermion mass eigenstates are superpositions of 
the couplings of the current eigenstates. 



a,f3=L,R 



(4.26) 



The elements of the 4x4 matrix T are given in Tab. 3a. The couplings dHkUf' ' 
current eigenstate basis fa,i3 = fL,R [normalized to 2(-\/2G'p')^/^] may be written as [|6|,p5| 





2 <I> 




2 <I> 











,2 



(4.27) 



for the neutral Higgs bosons, Hk = h,H,A. T3 = ±1/2 is the isospin of the [left- 
handed] sfermion and e/ its electric charge. The coefficients gf are given in Tab. 3b; in 
the decouphng limit, the coefficients gf reduce to 



cos(/5 + a) 
sin(/5 + a) 



sin 2/5 
-cos 2/? 



(4.28) 



[for the other coefficients, see eqs.(4.2)]. For the charged Higgs bosons, the couplings [also 
normalized to 2(-\/2G'i7')^/^] are 



9H+Ucda 



\/2 



gf + 



2 a0 



(4.29) 



with the coefficients (^Jy^ with a, P = L, R listed in Table 3c. 
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LL 


RR 


LR 


RL 


11 


008 6*/ COsOf/ 


sin 9f sin 9fi 


cos9f sin 6* J/ 


sin 9f cos9f' 


12 


— COS 9f sin 9f' 


sin 6'/ cos9f' 


cos 6*/ cos 9f' 


— sin 9f sin 9fi 


21 


— sin 9f cos 


cos ^/ sin 6*// 


— sin 9f sin 6*// 


cos9f cos9f/ 


22 


sin sin Of 


cos9f cos9f' 


— sin 9f cos9f/ 


— cos9f sin9f/ 



Tab. 3a: Transformation matrix for the Higgs couplings to sfermions in the presence of mixing. 
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9i 


^2 


^3 






h 


cos a/ sin [3 


— sin(a + /5) 


— sin a 1 sin /5 


cos a/ sin /3 


u 


H 


sin a 1 sin (5 


cos(a + /5) 


cos a/ sin /5 


sin a 1 sin /5 




A 








1 


-l/tg/9 




h 


— sin a/ cos /3 


— sin(a + /3) 


cos a/ cos /5 


— sin a/ cos /? 


d 


H 


cos a 1 cos /3 


cos(a + /3) 


sin a 1 cos /3 


cos a 1 cos /3 




A 








1 


-tg/9 



Tab. 3b: Coefficients in the couphngs of neutral Higgs bosons to sfermion pairs. 



LL 
9\/2 


9l/2 


LR 
9l/2 


RL 
9l/2 


m^/tan/3 + m^tan/3 
- sin 2(3 


rUumdiigP + l/tg/3) 



rudiid + AdtgP) 



rriuifi + Au/tgP) 




Tab. 3c: Coefficients in the couplings of charged Higgs bosons to sfermion pairs. 

Mixing between sfermions occurs only in the third-generation sector. For the first 
two generations the decay pattern is rather simple. In the limit of massless fermions, the 
pseudoscalar Higgs boson does not decay into sfermions since by virtue of CP-invariance 
it couples only to pairs of left- and right-handed sfermions with the coupling proportional 
to TUf. In the asymptotic regime, where the masses Mh^h± are large, the decay widths of 



the heavy CP-even and charged |^ Higgs bosons into sfermions are proportional to 



w 



sin^ 2(3 



(4.30) 



H 
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These decay modes can be significant only for low values tg/? [which implies sin^2/5 ~ 1]. 
However, in this regime the decay widths are inversely proportional to Mh, and thus 
cannot compete with the decay widths into charginos/neutralinos and ordinary fermions 
which increase with increasing Higgs mass. Therefore, the decays into first and second 
generations are unlikely to be important. 

In the case of the third generation squarks, the Higgs decay widths can be larger by 
more than an order of magnitude. For instance the decay widths of the heavy neutral 
Higgs boson into top squarks of the same helicity is proportional to 

T^H ^ ti) ^ (4.31) 

in the asymptotic region, and it will be enhanced by large coefficients [for small tg/3] 
compared to first/second generation squarks. Conversely, the decay widths into sbottom 

quarks can be very large for large tg/3. Furthermore, the decays of heavy neutral CP-even 
and CP-odd Higgs bosons to top squarks of different helicities will be proportional in the 
asymptotic region [and for the CP-even, up to the suppression by mixing angle] to 

r(if, A it) ~ [/X + A,/tan/3]^ (4.32) 

For II and At values of the order of the Higgs boson masses, these decay widths will 
be competitive with the chargino/ncutralino and standard fermion decays. Therefore, if 
kinematically allowed, these decay modes have to be taken into account. 

4.4 Numerical results 

The decay widths of the H, A and Higgs bosons into the sum of charginos and neu- 
tralinos, squark or slepton final states, as well as the standard and the total decay widths 
are shown in Figs. 6a, 7a and 8a as a function of mi/2 for two values of the pseudoscalar 
Higgs boson mass Ma — 300 and 600 GeV, and for positive and negative values; tg/3 is 
fixed to 1.75. 

Fig. 6a shows the various decay widths for the heavy CP-even Higgs boson. For Ma — 
300 GeV, the H ^ ti channel is still closed and the decay width into standard particles is 
rather small, being of (9(250) MeV. In this case, the decays into the lightest stop squarks 
which are kinematically allowed for small values of mi/2 will be by far the dominant decay 
channels. This occurs in most of the mi/2 range if > 0, but if < only for mi/2 ^ 50 
GeV which is already ruled out by CDF and LEP1.5 data. 

The decays into charginos and neutralinos, although one order of magnitude smaller 
than stop decays when allowed kinematically, are also very important. They exceed the 
standard decays in most of the mi/2 range, except for large values of mi/2 and /i < 
where no more phase space is available for the Higgs boson to decay into combinations 
of the heavy and light chargino/neutralino states. For small mi/2 values, chargino and 
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neutralino decays can be larger than the standard decays by up to an order of magnitude. 

As expected, the decay widths into sleptons are rather small and they never exceed the 
widths into standard particles, except for large values of mi/2- Note that due to the isospin 
and charge assignments, the coupling of the H boson to sneutrinos is approximately a 
factor of two larger than the coupling to the charged sleptons. Since all the sleptons of 
the three generations are approximately mass degenerate [the mixing in the f sector is 
very small for low values of tg/?], the small decay widths into sleptons are given by the 
approximate relation: T(H uu) ~ AV{H IlIl) — 4T{H — » IrIr). 

For larger values of Mh, Mh > 350 GeV, the decay widths into supersymmetric 
particles have practically the same size as discussed previously. However, since the H ^ tt 
channel opens up, the decay width into standard model particles becomes rather large, 
O{10 GeV), and the supersymmetric decays are no longer dominant. For Mh — 600 GeV, 
Fig. 6a, only the H ^ qq decay width can be larger than the decay width to standard 
particles; this occurs in the lower range of the mi/2 values. The chargino/neutralino 
decays have a branching ratio of ~ 20%, while the branching ratios of the decays into 
sleptons are below the 1% level. 

Fig. 6b and 6c show the individual decay widths of the heavy H boson with a mass 
Mh — 600 GeV into charginos, neutralinos, stop quarks and sleptons for the set of pa- 
rameters introduced previously. For decays into squarks, only the channels H — > iiii, iii2, 
and in a very small range of mi/2 values the channel i? — > are allowed kinematically 
[see Fig. 3c]. The decay into two different stop states is suppressed by the [small] mixing 
angle, and due to the larger phase space the decay H titi is always dominating. 

For the decays into chargino and neutralinos, the dominant channels are decays into 
mixtures of light and heavy neutralinos and charginos, in particular H — > X1X2 ^^'^ 
H — ^ X1X3 or X2X\- This can be qualitatively explained, up to phase space suppression 
factors, by recalling the approximate values of the relative branching ratios in the large 
1/^1 limit given in eqs.(4.18-20): T[H — > X1X2 ) ~ ^1 while T[H — >• X2X3) ~ 1 ^-^d 
T{H — > x?X3) ~ tdiV? 9w because sin 2/3 is close to one. The mixed decays involving X4 
are suppressed since they are proportional to (1 — sin 2/3), and all other decay channels 
are suppressed by powers of M|///^ for large values. 

The decay widths for the pseudoscalar Higgs boson are shown in Fig. 7a. There are no 
decays into sleptons, since the only decay allowed by GP-invariance, A fif2, is strongly 
suppressed by the very small f mixing angle. For M^ = 300 GeV, the decay into the 
two stop squark eigenstates, A tit2, is not allowed kinematically and the only possible 
supersymmetric decays are the decays into charginos and neutralinos. The sum of the 
decay widths into these states can be two orders of magnitude larger than the decay width 
into standard particles. 

For values of Ma above the threshold, the decay width into charginos and neutralinos 
is still of the same order as for low Ma, but because of the opening of the A ^ tt mode. 
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the total decay width increases dramatically and the chargino/neutralino decay branching 
ratio drops to the level of 20%. As in the case of the heavy CP-even Higgs boson H, the 
relative decay widths of the pseudoscalar boson into charginos and neutralinos, Fig. 7b, 
are larger for the channels involving mixtures of light and heavy neutralinos or charginos; 
the dominant decay modes are, roughly, A X1X2 ^-^^ A — > X1X4 X2X\- Again, this 
can be qualitatively explained, up to phase space suppression factors, by recalling the 
approximate formulae of eqs.(4.18-19), since the situation is the same as for H, with the 
two neutralino states x\ ^iid ^4 being interchanged. 

For small values of the common gaugino mass, mi/2 ^ 100 GeV, the decay mode of 
the pseudoscalar Higgs boson into stop squarks, A —>■ ^1^2, is phase space allowed. In 
this case, it is competitive with the top-antitop decay mode. As discussed previously, the 
1/M| suppression [and to a lesser extent the suppression due to the mixing angle] of the 
A —>■ tit2 decay width compared to r{A —>■ tt) will be compensated by the enhancement 
of the Atit2 coupling for large values of /i and At. 

Fig. 8a shows the decay widths for the charged Higgs boson. Since the dominant decay 
channel th is already open for Mij± ^ 300 GeV [although still slightly suppressed 

by phase space] , the charged Higgs decay width into standard particles is rather large and 
it increases only by a factor of ~ 4 when increasing the pseudoscalar mass to = 600 
GeV. The situation for the supersymmetric decays is quite similar for the two masses: the 
chargino/neutralinos decay modes have branching ratios of the order of a few ten percent, 
while the branching ratios for the decays into sleptons, when kinematically allowed, do 
not exceed the level of a few percent, as expected. Only the decay — > the only 
squark decay mode allowed by phase space [see Fig. 3c] for relatively low values of mi/2, 
is competitive with the tb decay mode. 

The decay widths of the charged Higgs into the various combinations of charginos and 
neutralinos are shown in Fig. 8b for Mjj± ~ 600 GeV. The dominant channels are again 
decays into mixtures of gauginos and higgsinos, since is large. The pattern follows 
approximately the rules of eq.(4.22), modulo phase suppression. 

As discussed in section 2, since the chargino, neutralino and sfermion masses scale as 
Ma, the situation for even larger values of the pseudoscalar Higgs boson mass. Ma ~ 1 
TeV, will be qualitatively similar to what has been discussed for Ma ~ 600 GeV. The only 
exception is that there will be slightly more phase space available for the supersymmetric 
decays to occur. 

5. Final Decay Products of the Higgs Bosons 

In this section, we will qualitatively describe the final decay products of the produced 
Higgs bosons. Assuming that Ma is large. Ma ^ 500 GeV, the decays into standard 
particles [and more precisely, the tt for the neutral and the tb decays for the charged Higgs 
bosons] always have substantial branching ratios, even for the value tg/3 = 1.75 which will 
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be chosen for the discussion. Therefore, to investigate decays into SUSY particles in the 
main production processes, e^e~ — > HA and H~^H~, one has to look for final states where 
one of the Higgs bosons decays into standard modes while the other Higgs boson decays 
into charginos, neutralinos or stop squarks. As discussed previously, the decays into the 
other squarks are disfavored by phase space, while the branching ratios into sleptons are 
always small and can be neglected. 

Let us first discuss the case where one of the Higgs bosons decays into chargino and 
neutralino pairs. 



e+e- ^ HA ^[u][x^x'] and [tt][x'x 



e e 



H^H- ^[tb][x^x'] (5.1) 



The lightest chargino Xi and next-to-lightest neutralino X2 decay into [possibly virtual] 
W, Z and the lightest Higgs boson h, assuming that decays into sleptons and squarks are 
kinematically disfavored. In the limit of large the decay widths [in the decoupling 



limit] are proportional to ^7 



^ixt^XiW^) ~ sin2 2/? (5.2) 
T{x',-xlZ) ~ cos2 2/?[(M2-Mi)/2/i]2 

Tixl-^Xih) ~ sin' 2/3 (5.3) 

In most of the parameter space, the W/Z/h are virtual [in addition to the three-body 
phase space factors, the decay widths are suppressed by powers of M2M^//i^] except 
near the upper values of mi/2. In the case of X25 the channel X2 ~^ Xi^ mode is always 
dominant although suppressed by additional powers of M| / /x^ compared to the X2 ~^ ^Xi 
mode, since both h and Z are off-shell, and the Z boson width is much larger than the 
width of the h boson for small values of tg/5. The radiative decay X2 ~^ Xil should play 
a marginal role except for very small values of mi/2 where the difference between the X2 
and Z boson masses becomes too large. 

For large values of mi/2, the sleptons become rather light compared to the gauginos and 
the decays of the light chargino and neutralino into leptons+sleptons are kinematically 
possible. In this case, these cascade decays become dominant since the partial widths for 
large are given by 

E nxl - II) = E 2r(xf - = '-^^M, (5.4) 
I I V2Tr 

and therefore not suppressed by powers of M^M2//i^, unlike the previous decay modes 
[we assume of course that there is no suppression by phase-space] . The sleptons will then 
decay into the LSP and massless leptons, leading to multi-lepton final states. 

The heavier chargino, in the absence of squark and slepton decay modes, will decay 
preferentially into the lightest chargino and neutralinos plus gauge or light Higgs bosons. 
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The decay widths, in units of CpM"^]!!]/ {S\/2tx) may be approximated in the decouphng 
hmit by ^ 

xt ^ xtz : r = i 

Xth : r = l 

^ xlW+ : r = 1 (5.5) 

The branching ratios for the various final states are roughly equal. Since xt is almost 
higgsino-like, the decay widths into sleptons and partners of the light quarks, when kine- 
matically allowed, are extremely small since they are suppressed by powers of m'j/M^. 
Because of the large value, only the decays into stop squarks and bottom quarks will 
be very important. This decay is allowed in most of the parameter space for M4 > 600 
GeV and, up to suppression by mixing angles, it is enhanced by a power 

r{xt W, Z, h) ^ sin^ /?(3 + tan2 Ow) ^ 

compared to the other decays. Therefore, when kinematically possible, this decay will be 
the dominant decay mode of the heavy charginos. 

For the heavier neutralinos, X34) the decay widths into W/Z/h bosons, again in units 
of GFM^|^|/(8y27r) may be be written in the decoupling limit as |3 

X3/4 ~ 





: r 


= ^tan2^vK(l ±sin2/5) 


X\h : 


r 


= itan^ ^14/(1 T sin 2/3) 


xlz 


: r 


= i(l±sin2/5) 


xlh : 


r 


= ^(lTsin2/3) 


xtw- 




r = 2 



(5.7) 

The dominant mode is the charged decay, X34 xt'^'y followed by the modes involving 
the h{Z) boson for X4(X3)- Because sin 2/3 ~ 1, only one of the h or Z decay channels is 
important. Here again, because of the higgsino nature of the two heavy neutralinos, the 
decay widths into sleptons and the scalar partners of the light quarks are negligible; the 
only important decays are the stop decays, X3 4 ^ ttiy when they are allowed kinematically 
[i.e. for not too large values of mi/2]. The ratio between stop and W/Z/h decay widths, 
up to suppression by mixing angles, is also given by eq.(5.6), and the stop decays will 
therefore dominate. 

We now turn to the case where one of the produced Higgs particles decays into stop 
squarks 

e+e" HA ^ [ti][iiii] and [ti][iit2] 
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e e 



[tb][hbi 



(5J 



From the squark mass plots, Fig. 3c, the only decay modes of the lightest stop squark 
allowed by phase space are 



txl , , bxt 



(5.9) 



Only the last decay mode occurs for relatively small values of mi/2, since m^^ < mt + m^a ^ 

in this case. For larger values of mi/2, h is heavy enough to decay into top quarks plus 
the lightest neutralinos. For these mi/2 values, the three decay modes of eq.(5.9) will have 
approximately the same magnitude since the chargino and the neutralinos are gaugino- 
like and there is no enhancement due to the top mass for the ti —>■ tx^ decays. 

The heavier stop squark, in addition to the previous modes, has decay channels with 
ti and Z/ h bosons in the final state 



^2 — tiZ , tih 



(5.10) 



These decays, in particular the decay into the lightest Higgs boson h, will be dominant 
in the large limit, since they will be enhanced by powers of /i^. 



Appendix A: Chargino and Neutralino Masses and Couplings 

In this Appendix we collect the analytical expressions of the chargino and neutralino 
masses and couplings, and we discuss the limit in which the Higgs-higgsino mass param- 
eter is large. 

The general chargino mass matrix [0, 



M 



c 



M2 V2Mw sin P 

y/2Mw COS (3 fi 



(Al) 



is diagonalized by two real matrices U and V, 

U = 0_ and V = 
where a is the matrix 



U*McV-^ 



0+ if detA^c > 
aO+ if detMc < 



a 



±1 
±1 



(A2) 



(A3) 



with the appropriate signs depending upon the values of M2, fi, and tan/5 in the chargino 
mass matrix. 0± is given by: 



O 



± 



cos 9± sin 9± 
— sin 9± cos 9± 



(A4) 
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with 



tan 26. 
tan 29, 



2V2Mw{M2 cos/5 + /i sin p) 

M| -/x2-2M^ cos /3 
2^/2Mw{M2 sin /3 + cos P) 



M|-/i2 + 2M2,cos/? 
This leads to the two chargino masses, the Xi.2 masses 

^xU= -^[mI + ^' + 2M^ 

T { {Ml - + 4M^ cos^ 2(3 + AMI,{MI + /i^ + 2M2/i sin 2(3) } 
In the hmit ^ M2, M^, the masses of the two charginos reduce to 



(A5) 



'(A6) 



A2 



Mo 



/i2 



(M2 + /isin2/?) 



/Vf2 

1/^1 + (M2sin2/5 + /i) 



(AT) 



where is for the sign of yU. For — ^ cx), the hghtest chargino corresponds to a pure wino 

M2, while the heavier chargino corresponds to a pure higgsino 

l/^l- 



state with mass ~ 
state with a mass m + = 



In the case of the neutralinos, the four- dimensional neutralino mass matrix depends on 
the same two mass parameters fj, and M2, if the GUT relation Mi = | tan^ 9w M2 — \M2 
[p!8[] is used. In the {—iB, —iW^, ifg) basis, it has the form 



M 



N 



Ml -Mzsw cos p Mzsw sin p 

M2 Mzcw cos P -Mzcw sin P 

-MzSw cos P MzCw cos P —fj, 

MzSw sin P —MzCw sin P —jj, 



(AJ 



It can be diagonalized analytically |^ by a single real matrix Z\ the [positive] masses 
of the neutralino states m^o are given by 



eim^o 



esm^o 



Ci + 



1 

r- 






1 

— a 
2 


1 

-a — 
2 




1/2 

1 - 


1 

-2" 


1 

-a — 
2 




1/2 

) - 


1 

— a 
2 


1 

—a — 
2 




1/2 

1 H- 


1 

— a 
2 



3 
1 

~ 3 



- 8^2/3)1/2 
C3 

^3a- 8^2/3)1/2 

(8a -8^2/3)1/2 
C3 



1/2 



.1/2 



nl/2 



1/2 



^C2/3)l/2 



(A9) 
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where Cj = ±1; the coefficients Ci and a are given by 

C, = {Mi + M2)/A 

M1M2 -Ml- - QCl 

2Ci [C2 + 2Cl + 2/^2] + M|(Mic^ + M2S^) - sin 2/5 
C1C3 - C^Ca - - MiMsAt^ + (Mic^^. + M^s^^) M^^ sin 2 p 



C2 
C3 
C4 



and 



2V3 



Re 



.2^xl/2]V3 
^ + ^27) 



(AlO) 



(All) 



with 



S 
D 



2 2 3 ^ 

C3 + ~ TtC'2C4 



27 



(C| + 12C4)^ - 27 



(A12) 



In the hmit of large values, the masses of the neutralino states simplify to 



m^o ~ Ml — 



m^o ~ M2 — 



Ml 
Ml 



(Ml + /xsin2/?) 



AS 



1M| 



(M2 + //sin2/3)c^ 



A*| + ^^£^(1 - sin 2/3) (// + MsS^ + Mic^) 



2 /i 

IMi 



m^o ~ |^| + -^e^(l + sin2/?)(^-M2S^-Mic^) 



(A13) 



Again, for — > 00, two neutralinos are pure gaugino states with masses m^o ~ Mi , 
m^o = M2, while the two others are pure higgsino states, with masses m^o ~ m^o ~ |//|. 

The matrix elements of the diagonalizing matrix, Z^j with i,j = 1, ..4, are given by 



Zil 

Zi4: 

Zil 



1 + 



H2 



Zis 

Zil J \Zii 
I Ml — eiUi^o 



+ 



Zi. 



iA 



Zil 



-1/2 



(A14) 



tan 9w M2 — eiUi^o 

Pl{Mi - eim^o){M2 - eiUi^o) - M| sin/3cos/3[(Mi - M2)c^ + M2 - e^m^o] 



Mz(M2 — eim^o)sw[ficos(3 + ejm^o sin (3) 
—eim^o{Mi — ejm^o)(M2 — ejm^o) — M| cos^/3[(Mi — M2)c^ + M2 — ejm^o] 



M2(M2 — eim^o)svK[A* cos/3 + eim^q sin /3) 

where is the sign of the ith eigenvalue of the neutralino mass matrix, which in the large 
|//| limit are: ei = £2 = 1 and £4 = —63 = e^. 
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Appendix B: Sfermion Masses and Mixing 



We now present the explicit expressions of the squark and slept on masses. Wc will assume 
a universal scalar mass mo and gaugino mass mi/2 at the GUT scale, and we will neglect 
the Yukawa couplings in the RGB's [see Appendix C]. For third generation squarks this 
is a poor approximation since these couplings can be large; these have been taken into 
account in the numerical analysis. 

By performing the RGE evolution to the electroweak scale, one obtains for the left- and 
right-handed sfermion masses at one-loop order [we include the full two-loop evolution 
of the masses in the numerical analysis] 



3 

H = ^0 + E W)ml/^ ± {Tsf - efsl,)Ml cos 2(3 

i=l 



(Bl) 



T^f and ej are the weak isospin and the electric charge of the corresponding fermion /, 
and Fi are the RGE coefficients for the three gauge couplings at the scale Q ~ M^, given 

by 



Ciif) 



(B2) 



The coefficients ftj, assuming that all the MSSM particle spectrum contributes to the 
evolution from Q to the GUT scale Mq, are 



6i=33/5 , 62 = 1 , &3 = -3 



(B3) 



The coefficients c(/) = (ci, C2, c^){f) depend on the hypercharge and color of the sfermions 
[Fl — Ll or Ql is the slepton or squark doublet] 



I 3/10 \ 




( 6/5 \ 


vn 










^ ; 



1/30 

c{Ql) =1 3/2 I , c{Ur) 
8/3 




ciDn) 




With the input gauge coupling constants at the scale of the Z boson mass 
ai{Mz) ~ 0.01 , a2{Mz) ~ 0.033 , a^iMz) ~ 0.118 
one obtains for the GUT scale Mq and for the coupling constant aa 

Mg ~ 1.9 X 10^^ GeV and aa = 0.041 
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(B4) 



(B5) 



(B6) 



Using these values, and including only gauge loops in the one-loop RGE's, one obtains 
for the left- and right-handed sfermion masses p9l 





= ml 


+ 6.28m2 ,„ 

1/ z 


+ 0.35M| 


cos(2/5) 




= ml 


+ 6.28m?/2 

_L/ Z 


- 0.42M| 


cos(2/5) 


"J? 


= ml 


+ 5.87m?/2 


+ 0.16M| 


cos(2/3) 




= ^1 


+ 5.82mi/2 


- 0.08M| 


cos(2/3) 




= ml 


+ 0.52m^/2 


+ 0.50M| 


cos(2/3) 




= ml 


+ 0.52m^/2 


- 0.27M| 


cos(2/?) 


^Ir 


= "^0 


+ 0.15mJ/2 


- 0.23M| 


cos(2/?) 



(B7) 

In the case of the third generation sparticles, left- and right-handed sfermions will 
mix; for a given sfermion f = t,b and f, the mass matrices which determine the mixing 
are 



ml + m 
where the sfermion masses m 



4 

JL 



m/(A/-/ir/) m^ +m| 



(B8) 



6 r 



are given above, m^ are the masses of the partner 
fermions and Vf, = = l/r^ = tg/5. These matrices are diagonalized by orthogonal 
matrices with mixing angles 6f defined by 

sin2^, = ^!^^#I^ , cos2^^, 



2 2 
ml — ml 

Jl 72 



'■7. 

2 2 
ml — ml 

Jl J2 



(B9) 



and the masses of the squark eigenstates given by 

1 



m 



/l,2 



m} + 



2 I 2 
m r + m r =F 



m 



m 



)2 + imj^Af - iirff 



(BIO) 



Appendix C: Renormalization Group Equations 

Finally, we collect for completeness the renormalization group equations for the soft- 
SUSY breaking parameters [the trilinear couplings, scalar masses as well as for /i and B] , 
including the dependence on At, A), and A^-. We restrict ourselves to the one-loop RGE's 
and we keep only the leading terms in the mass hierarchy in the MSSM with three fermion 



generations. The complete expressions for the RGE's can be found in Refs. [13,19 



For the trilinear couplings of the third generation sfermions, the RGE's are given by 



dt 
dAb 

~dr 

dAr 

~dr 



167r2 
2 



Y,c,g^^M, + Q\lAt + \lA, 



167r2 
2 



^Q^-^{Y.^^9^M^ + 3XlA, + AXlA, 



(CI) 
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while for the scalar masses of the third generation sfermions, one has 

^ - -^{-^9M-S9lMl-^-^glMi + X',X, + XlX,) 



dt 167r2 V 15^ 

dM^ 2 / 16 , , 16 , , , N 

dM^^ _ 2 / 4 16 



( - ^glM^, - ^glMi + 2A^X,) 



dt 167r2V 15^^ ' 3 

^ - ^(-|9?M?-3,M + A?X,) 
dMl 2 / 12 



( - ^9lM! + 2XlX^) 



dt 167r2 
The evolution parameter is defined hy t = \og{Q/MG), 

hi = ( 33/5 ,1,-3) 
Ci = (13/15, 3, 16/3) 
c[ = (7/15,3,16/3) 



(9/5,3,0) 



and 

X, = Ml^ + Ml + Ml^+Al 

X, = m1^+m4 + m|^+a^ 

For the first and second generation sfermions, these expressions reduce to 



and 



dAy, 
dt 


2 

167r2 


{^c,g^M, + X',A,) 




dAd 

dt 


2 

167r2 


[j24g^M, + XlA, + ^ 




dAe 
dt 


2 

167r2 






dMl 
dt 




- ^/^M! - 3gM - 




dt 


' f 

16^2 V 






dMj 
dt 


2 ( 

167r2 V 






dMl 
dt 


' ( 
167r2V 


- ^gfM^ - SglMi) 




dMl 
dt 


' ( 
167r2V 
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For the gauge coupling constants and the other soft-SUSY breaking parameters, the 
RGB's are given by 

dm|^^ _ 2/3 



( - -glM^, - ?,glMl + SA^X^ + A^X,) (Cll) 



dt 167r2 V 5 

dm|j2 _ 2/3 
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Fig. 1: Masses of the CP-even Higgs bosons h, H and of the charged Higgs particles H"^ 
as a function of Ma for two values of tg/3 = 1.75 and 50; the common squark mass Ms at 
the weak scale is fixed to Ms — 1 TeV and we take ji — At — Q. 
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Fig. 2: The correlation between mo and mi/2 for tg/9 = 1.75 and three values of — 
300, 600 and 900 GeV. The non-solid hnes show the boundaries which can be excluded by 
including the experimental bounds from LEP1.5 and Tevatron. 
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Fig. 3a: The masses of the Higgs bosons as a function of mi/2 for tg/3 = 1.75, for the 
two values uiq = 100 and 500 GeV and both signs of /i. 
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Fig. 3b: The masses of the two charginos (dashed hnes) and the four neutrahnos (sohd 
hnes) as a function of mi/2 for tg/3 = 1.75, Ma = 300 and 600 GeV and for both signs of 
II. The chargins/neutrahnos are ordered with increasing masses. 
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Fig. 3c: The masses of the two stop (soUd hnes), sbottom (dotted hnes) and first/second 
generation squark (dashed lines) eigenstates as a function of ra\ji for tg/3 = 1.75, Ma = 
300 and 600 GeV and for both signs of /x. 
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Fig. 3d: The masses of the charged sleptons (sohd and dotted hnes) and the sneutrino 
(dashed hnes) of the three generations as a function of mi/2 for tg/3 = 1.75, Ma = 300 
and 600 GeV and for both signs of n. 
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Fig. 4: The correlation between ttiq and mi/2 for tg/? ~ 50, jj < 0, and two values 
of Ma = 300 and 600 GeV. The boundary contours correspond to tachyonic solutions, 
m? < 0, M| < and M| < at the tree-level. 
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Fig. 5a: Cross sections for tlic pair production processes e^e HA and e+e H^H 
as a function of ^/s for tg/3 = 1.75 (solid lines) and tg/5 = 50 (dashed lines) and three 
values of Ma = 300, 600 and 900 GeV. 
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Fig. 5b: Cross sections for the production processes e~^e HZ, e^e — > hA and 
e^e" — > Huu as a function of ^/s for tg/3 = 1.75 and the values Ma = 300 and 600 GeV. 
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Fig. 6a: Decay widths (in GcV) of tlie lieavy CP-cvcn Higgs boson H into charginos and 
neutralinos (dotted lines), squarks (dashed hnes), sleptons (dash-dotted hnes), standard 
particles (dott-long-dashed hnes) and the total decay widths (sohd hnes) as a function 
of mi/2 for tg/3 = 1.75, Ma — 300 and 600 GeV and for both signs of ji. 
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Fig. 6b: Partial decay widtlis (in GeV) of tlie lieavy CP-even Higgs boson H into 
all combinations of chargino and neutralino pairs [ij = XiXj] ^ ^ function of mi/2 for 
tg/3 = 1.75, Ma = 600 GeV and for both signs of pi. 
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Fig. 6c: Partial decay widtlis (in GeV) of the heavy CP-even Higgs boson H into stop 
and sbottom squarks and into slepton pairs as a function of mi/2 for tg/3 = 1.75, — 600 
GeV and for both signs of ji. 
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Fig. 7a: Decay widths (in GcV) of the pseudoscalar Higgs boson A into charginos and 
neutrahnos (dotted hnes), stop squarks (dashed hnes), standard particles (dott-long- 
dashed hnes) and the total decay widths (solid hnes) as a function of mi/2 for tg/3 = 1.75, 
Ma — 300 and 600 GeV and for both signs of 
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Fig. 7b: Partial decay widtlis (in GeV) of tiic pscudoscalar Higgs boson A into all 
combinations of chargino and neutralino pairs [ij = XiXj] ^ function of mi/2 for tg/3 = 
1.75, Ma = 600 GeV and for both signs of 
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Fig. 8a: Decay widths (in GeV) of tlie cliarged Higgs bosons into charginos and neutrali- 
nos (dotted lines), squarks (dashed hnes), sleptons (dash-dotted hnes), standard particles 
(dott-long-dashed lines) and the total decay widths (sohd hnes) as a function of mi/2 for 
tg/3 = 1.75, Ma = 300 and 600 GeV and for both signs of 
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Fig. 8b: Partial decay widtlis (in GeV) of tiic charged Higgs boson H"^ into all combi- 
nations of charginos and neutralinos [ij = Xi^Xj] ^ function of mi/2 for tg/3 = 1.75, 
Ma — 600 GeV and for both signs of 
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Production of Heavy Neutral MSSM Higgs Bosons 
a complete 1-loop calculation. 



V. Driesen, W. Hollik and J. RosiekQ 



Institut fiir Theoretische Physik, Universitat Karlsruhe, D-76128 Karlsruhe, Germany 



The complete 1-loop diagrammatic calculations of the cross sections for the 
neutral Higgs production processes e~^e~ — > Z^H^{Z^h^) and e'^e~ A^H^{A^h^) 
in the MSSM are presented and compared the with the corresponding results of the 
simpler and compact effective potential approximation. 

1. Introduction 

In order to experimentally detect possible signals of the neutral MSSM Higgs bosons, 
detailed studies for the decay and production processes of Higgs boson are required. As 
has been discovered several years ago [1-3], radiative corrections in the MSSM Higgs sector 
are large and have to be taken into account for phenomenological studies. Three main 
approaches have been developed to calculate the 1-loop radiative corrections to the MSSM 
Higgs boson masses, production and decay rates: 

a) The Effective Potential Approach (EPA) 0. 

b) The Renormalization Group approach (RGE) 

c) The diagrammatic calculation in the on-shell renormalization scheme (Feynman 
Diagram Calculation, FDC) The masses are calculated from the pole positions 
of the Higgs propagators, and the cross sections are obtained from the full set of 
1-loop diagrams contributing to the amplitudes, including 0]: 

- the most general form of the MSSM lagrangian with soft breaking terms, 

- the virtual contributions from all the particles of the MSSM spectrum, 

- all 2-, 3- and 4-point Green's functions for a given process with Higgs particles, 

- the momentum dependence of the Green's functions, 

- the leading reducible diagrams of higher orders corrections. 

The experimental searches for Higgs bosons at LEPl and studies for the fu- 
ture searches at higher energies [j^ conventionally make use of the most compact ef- 
fective potential approximation. We present the complete 1-loop diagrammatic results 



compact EPA approximation and discuss the typical size of the differences. 



* Supported in part by the Alexander von Humboldt Stiftung and by the Polish Committee for Scientific 
Research. 
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2. Outline of the calculations 

The tree level potential for the neutral MSSM Higgs bosons can be written as: 

= mlHl + mlHl + e.,{m\2H{Hi + H.c.) + ^^^^{Hl - H^f + ^^{H,H,f (1) 

Diagonalization of the mass matrices following from the potential (|l|) leads to three phys- 
ical particles: two CP-even Higgs bosons , and one CP-odd Higgs boson , and 
defines their tree-level masses itih , Tn^ and , with itlh > Tn^ , and the mixing angles 
a, (5. The way of calculating the radiative corrections in the EPA and FDC methods is 
briefly described as follows: 

In the EPA, the tree level potential V^^^ is improved by adding the 1-loop terms 0: 

V")W=)=V'<">W^) + ji5 Y. StrM'Ls^^l) (2) 

quarks \ * / 

squarks 



Figure 1: Classes of diagrams contributing to the e+e — > Z'^hP[H^) process in the FDC 
approach. 

where V'^'^\Q'^) is the tree level potential evaluated with couplings renormalized at the 
scale Q^, and Str denotes the supertrace over the third generation of quark and squark 
flelds contributing to the generalized mass matrix AA^. The 1-loop potential V^*-^^ is 
rediagonalized yielding the 1-loop corrected physical masses Mh , and the effective 
mixing angle cte// (for explicit formulae see 0). 

In the FDC the 1-loop physical Higgs boson masses are obtained as the pole positions of 
the dressed scalar propagators. Mfj and M| are given by the solution of the equation (Q). 
For the calculations of the cross sections we need the full set of 2-, 3- and 4-point functions. 
In Fig. |l| the diagrams contributing to the e^e^ — ^ Z^h^{H^) process are collected. The 
diagrams contributing to the e^e~ — * A^h^{H^) process can be obtained by changing 
into on the external line and skipping the diagrams i), j). 



Re 



[p' - ml - {p' -ml- J^HHip')) - SL(P')] = (3) 



The formulae for the cross sections obtained in the FDC differ from the Born expres- 
sions, because not only the effective masses are corrected but also new form factors and 
momentum dependent effects are considered (see M for a detailed description). 



80 



3. Results on production cross sections 



In this section we present tiie results for Z^H^{Z^h^) and A^H^{A^h^) production 
from the FDC and discuss the quality of simpler EPA approximation. In all figures we 
use as an example the set of parameters listed in Table |l]. fi is the parameter describing the 
Higgs doublet mixing in the MSSM superpotential. M2 denotes the SU(2) gaugino mass 
parameter. For the U(l) gaugino mass we use the value Mi = | tan^ 6'vi/^2, suggested 
by GUT constraints. Mgq, Mgi, At and Ab are the parameters entering the sfermion mass 
matrices (for the detailed expressions see e.g. [^). For simplicity we assume a common 
value Msq for all generations of squarks, and a common Mgi for sleptons. 



Parameter 




Ma 


M,g 


M,i 


M2 




At = Ab 


Value (GeV) 


175 


200 


1000 


300 


1000 


500 


1000 



Table 1: Parameters used for the numerical analysis. 



From the theoretical point of view, the most convenient parameters for the Higgs 
sector are the mass M^ of the CP-odd Higgs boson and the ratio tan/5 = ^. From the 
experimental point of view it is more natural to use, depending on the process considered, 
the masses Mh or Mh of the CP-even Higgs instead of the formal quantity tan (3. 

As a first step, we apply the conventional Ma , tan/3 parameterization. Figs. ^ 
and ^ show the production cross sections for the processes a{e'^e~ — > Z^h^, A^h^) and 
a{e~^e~ Z^H^, A^H^) for = 500 GeV. For the chosen set of parameters the numer- 
ical differences can reach 30% at = 500 GeV. They become more important with 



Figure 2: Comparison of the cross sections a{e^e — > Z^h^,A^h^) obtained in the EPA and 
FDC. Parameters as given in Table 0, a/s = 500 GeV. 
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Figure 3: Comparison of the cross sections a{e~^e~ Z^H^,A^H^) obtained in the EPA and 
FDC. Parameters as given in Table 0, = 500 GeV. 

increasing energies, exceeding 40% at 1 TeV. Note, however, that in the region of large 
cross sections the EPA accuracy is better (20% at 500 GeV). More detailed discussion of 
the lighter CP-even Higgs boson production can be found in ref. 0. 

Fig. ^ shows the production cross sections for the processes a{e^e^ Z^H^, A^H^) 
as a function of . The effect of the additional form factors included in the FDC grows 
when center-of-mass energy increases. For ^Js = 1.5 TeV the differences between FDC 
and EPA can reach 50% for the a{e^e^ Z^H^) production channel. In addition, the 
angular dependence of the cross section given by the FDC is modified in compare to the 
effective Born approximation. 

We now turn to the more physical parameterization of the cross sections in terms of the 
two Higgs boson masses and or Mh ■ This parameterization is more clumsy in the 
calculations, but it has the advantage of physically well defined input quantities avoiding 
possible confusions from different renormalization schemes. Varying {Ma and other 
input quantities fixed) we obtain tan/? and ctzh, o'ah as functions of Mh ■ For the 
parameter values given in Table |I], the differences between the tan/5 values obtained in 
the EPA and FDC can reach 10% (up to 20% for smaller Ma ~ 100 GeV). Also significant 
differences can occur for the cross sections, as displayed in Fig. ^ where the predictions of 
EPA and FDC for the azH and a ah are plotted as functions of Mh . The typical size of 
differences between the methods is 10-20% for = 500 GeV, but they may became as 
large as 60% in case of the process a{e~^e~ — ^ Z^H^). This particularly large deviation 
occurs for large Mh values, corresponding to small tan/? < 1 (compare Fig. |). 

We have analyzed also the dependence of the differences between the EPA and the 
FDC predictions on the SUSY parameters: sfermion and gaugino masses, parameter 
and sfermion mixing parameters. In most cases the variation of those parameters does 
not have a large effect on the size of the differences between the EPA and FDC (a more 
detailed discussion can be found in ref. JTUf)- 

To give a more global impression of the typical size of the differences between the EPA 
and FDC results, we have chosen 1000 random points (for each ^Js value in Table |^) from 
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the hypercube in the MSSM parameters space with the following bounds: 



0.5 < tan/? < 50 

-500 GeV < ;U < 500 GeV 



50 GeV <Ma < 250 GeV 
200 GeV < Ms < 1000 GeV 



200 GeV < Msq = 2Msi < 1000 GeV -M,g < At = Ab < Msg 



We define the relative differences for the masses and cross sections as follows: 



6X 



EPA 



XFDC 



(4) 



where X can be chosen as , Mh , cfzh^ (^zh-, o'ah ot aAh- 

We calculated the quantities 5Mf (5M|^^, 5af^^, 5cr|f^, 5a|^^and ^af^^and 
averaged them (and also their absolute values) arithmetically over all generated points of 
the parameter space. The average mass differences are small and equal \5M^^^\ = 2% 
and \SM§^^\ = 1%. The results for the cross sections are summarized in Table ^ It 
shows that the predictions of both methods deviate in particular for azH- 



Vs (GeV) 




^^Ah 




^^AH 










500 


16.4% 


-2.4% 


57% 


4.4% 


21% 


31% 


62% 


14% 


1000 


10.3% 


1.1% 


56% 


-3.0% 


15% 


31% 


62% 


14% 


1500 


4.2% 


4.9% 


53% 


-9.0% 


17% 


32% 


63% 


18% 



Table 2: Differences between the EPA and FDC predictions averaged over a random sample of 
parameters. 

Summarizing, comparisons between the FDC predictions with the simpler EPA ap- 
proximation have shown that at = 500 GeV the EPA has an accuracy of typically 
10-20% in the parameter regions where the cross sections are large. The differences be- 
come larger with increasing energy, where also modifications of the Born-like angular 
distributions are more visible. The use of the physical input variables M^, Mh or M^, 
Mh avoids ambiguities from the definition of tan (3 in higher order, but the observed dif- 
ferences remain of the same size. For a better accuracy, the full FDC would be required. 

Recently some papers on the leading 2-loop corrections to the CP-even MSSM Higgs 
boson masses have been published [111]. The main conclusion is that 2-loop corrections 



are also significant and tend to compensate partially the effects of 1-loop corrections. The 
calculations are based on the EPA and RG methods. Since the main emphasis of this 
study is to figure out the difference between complete and approximate results in a given 
order, we have not implemented the 2-loop terms. They would improve the 1-loop FDC 
results in the same way as the approximations and thus do not influence the remaining 
differences which can only be obtained by an explicit diagrammatic calculation. 

The library of FORTRAN codes for the calculation of the 1-loop radiative corrections 
in the on-shell renormalization scheme to the MSSM neutral Higgs production and decay 
rates ^ can be found at the URL address: 

http:/ /itpaxpl . physik. uni-karlsruhe.de/r^rosiek/neutraLhiggs. html 
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Figure 4: Comparison of the cross sections a{e'^e Z^H^,A^H^) as a function 
^/s obtained in the EPA and FDC. tan/3 = 2, other parameters as given in Table |l|. 



Figure 5: Comparison of the cross sections a{e^e — > Z^H^,A^H^) as a function of Mh 
the EPA and FDC. Parameters as given in Table |l], ^/s = 500 GeV. 
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Abstract 

We discuss the one-loop electroweak corrections to the pair production of charged 
Higgs bosons e'^e" H^H~ in the Minimal Supersymmetric Standard Model. 

In contrast to hadronic machines, a high energy e+e~ collider in the TeV range will 
be a rather unique place to discover and study charged higgses in a clean environment. 
These would be produced either in pairs , our main concern here, or in associate (rare) 
production with W^. It was first found in that loop corrections from matter fermions 
and their susy partners (mainly the (t,fe), {h,t) sector), are likely to change the tree-level 
result at a/s = 500 GeV by asmuch as 10% dip in the cross-section. The effect could 
even lie between —25% and 25% and perhaps grow out of perturbative control, though 
in a reasonable range of the model-parameters. Such a sensitivity to loop effects appears 
to be related to the fact that at tree-level the 7 and Z mediated process is exclusively 
controlled by UsiX) x UwziX) gauge invariance and thus knows nothing about the non- 
standard extension whatsoever. 

The aim of the present study is to improve on the previous one by including: a) the 
complete Higgs sector contributions (self-energies, vertices and boxes), b) the infrared 
part, including initial and final soft photon radiation as well as 77 and 7Z boxes, c) The 
complete set of charginos/neutralinos/e/z/ box diagrams, and thus to identify the various 
origins of large effects, whether in the MSSM or in a type II two-Higgs-doublet model 
(THDM-II). 

It turns out that besides the sensitivity to the heavy quark-squark sector there are, on 
one hand large effects from the soft photon radiation and on the other, possibly important 
effects in the purely Higgs sector. The latter case occurs when deviations from the tree- 
level supersymmetric H^H~ —H^{h?) couplings are allowed bringing in increasingly large 
effects for increasing values of tan /? at a given y/s. To quantify such effects we thus allow 
for a general deviation from the supersymmetric relations among the bare parameters of 
the Higgs potential as follows : 

Ai = A2 + 5i2 , A3 = ^(/ + ^")-Ai + (53i 

o 

A4 = 2Ai - ^ 541 , A5 = -i(/ + ^?'=^)+2Ai + (55i (1) 
A6 = -^(/ + ^7")+2Ai + 56i 
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The X[s are as defined in and the (softly broken) susy case corresponds to 6i = 0. Eq.(l) 
translates in a definite way into deviations from the MSSM tree-level relations among the 
higgs masses, tan 2a and tan/5, as well as the couplings of the higgs sector, leading to 
6 free parameters. However, one can find conditions relating the 6's in such a way to 
preserve these relations even in a non-susy case. These conditions which we dub "quasi- 
susy" , are certainly not generic but constitute a good ground to test minimal deviations 
from the MSSM in a simple way, since one then has just one extra free parameter (ex. 
A3) besides tan (3 and Mh± in the higgs sector. In quasi-susy the only deviations from the 
MSSM tree-level triple Higgs couplings that contribute to one-loop order in our case are 
in if+ — H~ — {H^, h^). For large tan/3 these couplings behave as: 

{H+H-{H\ h''))susy - igMw{cos{(3 - a)A(i,2) + (1 T cos(2a))cos(«)tan(/?)A3 (2) 

the A's being functions of the 5's of eq.(l), and vanish in the MSSM. 

Fig.l illustrates how the Higgs sector contributions can counterbalance those of the 
heavy quarks found in for large tan(/5), but only near threshold. Far from thresh- 
old most of the effects become again negative, except for WW boxes. Furthermore the 
"neutral" model-independent contributions, including soft bremsstrahlung, obtained by 
adding one photon (or Z) line to the tree diagrams depend loosely on Mh± or ^/s and 
contribute at the level of —17% for Ai?^ ~ O.lE'feeam- In Fig. 2 we show (excluding those 
"neutral" contributions) the integrated cross-section for two values of Mh± and tan(/3). 
In THDM-II the total loop effect increases (negatively) with increasing tan(/?), the far- 
ther one goes from production threshold. In the MSSM (A3 = \zsusy ) the leading effects 
come exclusively from the heavy quark-squark sector and the conclusions of p[ remain 
unaltered in this case. [For instance the 150 boxes involving j')^ jclv largely cancel 
among each other leading at most to 1 — 3% effect for a wide range of sparticle masses.] 
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Fig.l: Contributions in % to the integrated cross section in quasi-susy, A3 = — 0.61,M^± = 
220i GeV; a) Higgs sector, tan/3 = 2; b) Higgs sector, tan/3 = 30; c) virtual Z, 7 and soft 
bremsstrahlung; d) virtual W boxes; e) matter fermion sector, mtop = 180Gey,tan/3 = 30; f) 
same as e) but with tan/3 = 2. 



Fig.2: a) Tree-level, Mh± = 220GeV; b) quasi-susy, A3 = -0.61 (MSSM value -0.71), tan/3 = 
30; c) quasi-susy, A3 = — 0.61, tan /3 = 2; d) Tree-level, Mfj± = 430GeV; e) quasi-susy, A3 = —2.6 
(MSSM value -2.84), tan/3 = 30; f) same as e) but with tan/3 = 2; rritop = ISOGeV. 
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Abstract 

The cross sections for the multiple production of the lightest neutral Higgs 
boson at high-energy e'^e~ colliders are presented in the framework of the 
Minimal Supersymmetric extension of the Standard Model (AiSSAi). We 
consider production through Higgs-strahlung, associated production of the 
scalar and the pseudoscalar bosons, and the fusion mechanisms for which 
we use the effective longitudinal vector-boson approximation. These cross 
sections allow one to determine trilinear Higgs couplings Xnhh and Xhhh, which 
are theoretically determined by the Higgs potential. 

1. Introduction 

The only unknown parameter in the Standard Model {SA4) is the quartic coupling of the 
Higgs field in the potential, which determines the value of the Higgs mass. If the Higgs 
mass is known, the potential is uniquely fixed. Since the form of the Higgs potential is 
crucial for the mechanism of spontaneous symmetry breaking, i.e. for the Higgs mech- 
anism per se, it will be very important to measure the coefficients in the potential once 
Higgs particles have been discovered. 

If the mass of the scalar particle is less than about 150 GeV, it very likely belongs 
to the quintet of Higgs bosons, h, H, A, predicted in the two-doublet Higgs sector 
of supersymmetric theories []I| [h and H are the light and heavy CP-even Higgs bosons, 
A is the CP-odd (pseudoscalar) Higgs boson, and H"^ is the charged Higgs pair]. The 
potential of the two doublet Higgs fields, even in the Minimal Supersymmetric Standard 
Model {M.SS}A), is much more involved than in the Standard Model 0. If CP is 
conserved by the potential, the most general two-doublet model contains three mass 
parameters and seven real self-couplings. In the M.SSM., the potential automatically 

* Supported by Deutsche Forschungsgemeinschaft DFG (Bonn). 
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conserves CP; in addition, supersymmetry fixes all the Higgs self-couplings in terms of 
gauge couplings. The remaining three free mass parameters can be traded in for the two 
vacuum expectation values (VEV's) of the neutral Higgs fields and one of the physical 
Higgs masses. The sum of the squares of the VEV's is fixed by the W mass, while the 
ratio of VEV's is a free parameter of the model called tan (3 . It is theoretically convenient 
to choose the free parameters of the J^SSA4 Higgs sector to be tan j3 and Ma, the mass 
of the CP-odd Higgs boson A. The other Higgs masses and the mixing angle a of the 
CP-even neutral sector are then determined. Moreover, since all coefficients in the Higgs 
potential are also determined, the trilinear and quartic self-couplings of the physical Higgs 
particles can be predicted theoretically. By measuring these couplings, the Higgs potential 
can be reconstructed - an experimental prima facie task to establish the Higgs mechanism 
as the basic mechanism for generating the masses of the fundamental particles. 

The endeavor of measuring all Higgs self-couplings in the AiSSA4 is a daunting task. 
We will therefore discuss a first step by analyzing theoretically the production of two light 
Higgs particles of the J^SSAi. These processes may be studied at the proton collider 
LHC [P] and at a high-energy e'^e~ linear collider. In this paper we will focus on the 
e'^e" accelerators that are expected to operate in the first phase at an energy of 500 GeV 
with a luminosity of about / £ = 20 fb~^, and in a second phase at an energy of about 
1.5 TeV with a luminosity of order J C = 200 fb~^ per annum [Q. They will allow us to 
eventually study the couplings XHhh and Xhhh- The measurement of the coupling Xhaa 
will be very difficult. 

Multiple hght Higgs bosons h can [in principle] be generated in the J^SSJ^ by four 
mechanisms^: 

(i) Decay of the heavy CP-even neutral Higgs boson, produced either by if-strahlung 
and associated AH pair production, or in the WW fusion mechanisms. Fig. la. 



e+e- ZH, AH 



H ^hh (1) 



Associated production e'^e — > hA followed by yl — > hZ decays leads to hhZ background 
final states. 

(ii) Double Higgs-strahlung in the continuum, with a final state Z boson. Fig. lb, 

e+e" ^ Z* ^ hhZ (2) 

(iii) Associated production with the pseudoscalar A in the continuum. Fig. Ic, 

e+e" ^ Z* ^ hhA (3) 

(iv) Non-resonant WW{ZZ) fusion in the continuum. Fig. Id, 

e+e" iy^iy^W*W* v^vM (4) 



^The production of two light Higgs bosons, e"'"e — > hh, through loop diagrams does not involve any 
trilinear Higgs coupling; the production rates are rather small 
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The cross sections for ZZ fusion in (1) and (4) are suppressed by an order of magnitude. 
The largest cross sections can be anticipated for the processes (1), where heavy on-shell 
H Higgs bosons decay into pairs of the hght Higgs bosons. [Cross sections of similar 
size are expected for the backgrounds involving the pseudoscalar Higgs bosons.] We have 
derived the cross sections for the four processes analytically; the fusion process has been 
treated in the equivalent particle approximation for longitudinal vector bosons. 

We will carry out the analysis in the M.SSM. for the value tan/3 = 1.5. [A summary 
will be given in the last section for all values of tan/3 ]. In the present exploratory 
study, squark mixing will be neglected, i.e. the supersymmetric Higgs mass parameter 
[I and the parameter At in the soft symmetry breaking interaction will be set to zero, 
and the radiative corrections will be included in the leading m\ one-loop approximation 
parameterized by 

or^ / ]\/r2\ 

(5) 




with the common squark mass fixed to Ms = 1 TeV. In terms of tan/3 and M^, the 
trilinear Higgs couplings relevant for our analysis are given in this approximation by 

6 COS (y 

Khh = 3cos2asin(/? + a) + 3 . (6) 

M| sm p 

^Hhh = 2sin2asin(/3 + a) — cos2acos(/? + a) + 3— -2"^: — ^ cos^ a 



■ M| sin p 



In addition, the coupling 



AftAA = COS 2/3 sin(/3 + a) + —j cos^ p (7) 

M| sm p 

will be needed even though it turned out - a posteriori - that it cannot be measured using 
the experimental methods discussed in this noteQ. As usual, these couplings are defined 
in units of (2\/2G'i?)^/^M|; the h, H, masses and the mixing angle a can be expressed 
in terms of Ma and tan/3 [see e.g. Ref. [§] for a recent discussion]. 

In the decoupling limit for large A, H and if^ masses, the lightest Higgs particle 
becomes iSA^-like and the trilinear hhh coupling approaches the SM. value Xhhh 
Ml/M^. In this limit, only the first three diagrams of Fig. lb and Id contribute and the 
cross-sections for the processes e^e~ — > hhZ and WW — hh approach the corresponding 



cross sections of the SAi WWiTl . 



2. H Production and hh Decays 

If kinematically allowed, the most copious source of multiple h final states are cascade 
decays H hh, with H produced either by Higgs-strahlung or associated pair production 

^For small masses the decay h — > AA could have provided an experimental opportunity to measure 
this coupling. However, for tan/3 > 1, this area of the M.SSM. parameter space is excluded by LEP 0|. 
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a[e e 



a[e e 



ZH) 
AH) 



967rs 



a] 



Xf[\z + l2Ml/s] 
(1 - Ml/sf 



X 



3/2 



(8) 
(9) 



The Z couplings to electrons are given by = — l,fe = — 1 + 4sin^6'vi/ and is the 
usual two-body phase space function Xj = (1 — Mf/s — Mfj/sY —AMjMfj/s"^. The cross 
sections (8) and (9) are shown in Fig. 2 for the total e~^e~ energies y/s = 500 GeV and 
1.5 TeV as a function of the Higgs mass Mh for a small value of tan j3 = 1.5 where the H 
cascade decays are significant over a large mass range. As a consequence of the decoupling 
theorem, associated AH production is dominant for large Higgs masses. 

The trilinear Hhh coupling can be measured in the decay process H ^ hh 



r{H hh) 



1 - 



H 



AMI 
Ml 



1/2 



(10) 



if the branching ratio is neither too small nor too close to unity. This is indeed the case, 
as shown in Fig. 3a, for H masses between 180 and 350 GeV and small to moderate tan (3 
values. The other important decay modes are WW* / ZZ* decays. Since the H couplings 
to the gauge bosons can be measured through the production cross sections of the fusion 
and Higgs-strahlung processes, the branching ratio BR(/7 hh) can be exploited to 
measure the coupling Xuhh- 

The ZH final state gives rise to resonant two-Higgs [hh] final states. The AH final 
state typically yields three Higgs h[hh] final states since the channel A —>■ hZ is the 
dominant decay mode in most of the mass range we consider. This is shown in Fig. 3b 
where the branching ratios of the pseudoscalar A are displayed for tan/3 = 1.5. 

Another type of two-Higgs hh final states is generated in the chain e~^e~ Ah —>■ 
[Zh]h, which does not involve any of the Higgs self-couplings. However, in this case, the 
two h bosons do not resonate while [Zh] does, so that the topology of these background 
events is very different from the signal events. The size of the e'^e" — >■ hA background 
cross section is shown in Fig. 2 together with the signal cross sections; for sufficiently 
large M4, it becomes small, in line with the decoupling theorem 

A second large signal cross section is provided by the WW fusion mechanism. [Since 
the NC couplings are smaller compared to the CC couplings, the cross section for the ZZ 
fusion processes in (1) and (4) is ~ 16cos^6'vi/, i-G. one order of magnitude smaller than 
for WW fusion.] In the effective longitudinal W approximation one obtains 



cr e^e 



G^M^ 



4V2 



TT 



1 + 



log 



Ml 



2 1 



Ml' 



cos2(/5-a) (11) 
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The magnitude of the cross section^ e~^e^ — > Huei^e is also shown in Fig. 2 for the two 
energies ^/s = 500 GeV and 1.5 TeV as a function of the Higgs mass Mh and for tan/5 = 
1.5. The signals in e^e~ [hh] + missing energy are very clear, competing only with 
iJ-strahlung and subsequent neutrino decays of the Z boson. Since the lightest Higgs 
boson will decay mainly into bb pairs, the final states will predominantly include four and 
six b quarks. 

At a/s = 500 GeV, about 500 signal events are predicted in the mass range of ~ 
200 GeV for an integrated luminosity of / £ = 20 fb~^ per annum; and at y/s = 1.5 TeV, 
about 8,000 to 1,000 signal events for the prospective integrated luminosity of / £ = 200 
fb~^ per annum in the interesting mass range between 180 and 350 GeV. Note that for 
both energies, the Ah background cross section is significantly smaller. 



3. Non-Resonant Double hh Production 

The double Higgs-strahlung e"'"e~ Zhh, the triple Higgs production process e"^e^ 
Ahh and the WW fusion mechanism e'^e^ Uei^ehh outside the resonant H hh 
range are disfavored by an additional power of the electroweak coupling compared to the 
resonance processes. Nevertheless, these processes must be analyzed carefully in order to 
measure the value of the hhh coupling. 



3.1 e+e- Zhh 

The double differential cross section of the process e~^e~ 
da{e~^e~ hhZ) 



dxidx2 



384^2- 



hhZ, Fig. lb, is given by 
A 



12) 



The couplings have been defined in the previous section. xi^2 = 2-Ei,2/a/^ are the scaled 
energies of the Higgs particles, X3 = 2 — xi — X2 is the scaled energy of the Z boson; 
Uk = I — Xk- The scaled masses squared are denoted by fii = Mf /s. In terms of these 
variables, the coefficient A in the cross section may be written as: 



A 



^ ^ sm^{(3 — a) 



f: 



f2 



X 



— + - 

yi+ fJ'h- IJ'A 1/2 + fJ'h- IJ'A^ 



yi + fj'h- fj'Z 

a 



y2 + jj'h- jj'Z 



+ 



cos''(/? — a) 



+ 



8fiz{yi + f^h - fiz) 



l^z 



+ 



sin^iP - a)U ^ cos^(/? - a)/6 



yi + 
h 



yi + i^h- /iz y2 + i^h- /iz 



yi 
{yi 



* 1/2} 



with 



a 



sm{(3 - a)Xhhh ^ cos{(3 - a)\Hhh 



y3 + fJ'Z - IJ'h 



Z/3 + Atz - fJ-H 



^ sin^(/?-a) ^ sin^(/3 - a) 



1 



(13) 



(14) 



yi + f^h- f^z y2 + f^h- f^z '^fj'Z 

■^In the effective W approximation, the cross section may be overestimated by as much as a factor of 2 
for smah masses and/or smah cm. energies. Therefore we display the exact cross sections [fsl in Fig. 2. 
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[omitting the small decay widths of the Higgs bosons]. Only the coefficient a includes 
the Higgs self-couplings Xuhh and \hhh- Introducing the notation j/o = — l/2)/2, the 
coefficients fi which do not involve any Higgs couplings, are defined by 



/o 
/i 

h 
h 
h 
h 

f7 

fs 



i + Vif -Ajiziyi + iJ'h) 



(15) 



yl 



2/iz(l/i + 



[2/iz(Aiz - 2/i/i + 1) - (1 + + ?/2)] [Aiz(/iz -y\-y2- + 2) - y^y-^ 
yl + - Z/i - 2/2 + yUz - 4/i/,) 1 + 2/1 + 2/2 + 2/o + - - 2?/i 



2/0 + - 2/1 - 2/2 + A^z - 4/i/, 



2/0^ 



2/i| 



2/1^ 



4/i^(|/l+2/i;,)+/i^[(l + |/l)(3l/i 
2 



1 + ^J'z{^^z - 2/1 - 2/2 - 4/2/, + 2) 
-l/2) + 2] -2/?(l + 2/i+2/2) 



2/12/2 

z - /^z(2/2 + 32/1 + 8/i,i - 2) + 2/iz2/o (1 + 2/i + 2/o) + 22/i2/o - 2/o(2/i + 2/2-2) 
[/iz(4/i/, - - 1 + 22/1 - 2/o) - 2/i2/o] [Aiz(4/i/i - - 1 + 32/i) - (1 + 2/o)(l + 2/i)] 
[/iz(4Aifc - /iz - 1 + 2?/i - yo) - 2/i2/o] [Aiz(4/i/i - - 2 + yi) + (1 - ?/o)(l + 2/i)] 



In the decoupling limit, the cross section is reduced to the SM. cross section for which 



A=-h+ 1 

2 ° 4/x|(?/i + /i/, - yUz) 

with the /j's as given above, and 



/: 



/2 



2/1 + /W/i - A^z 2/2 + Ai/i - /iz 



+ {2/1 ^ 2/2} 



1 

a = -- 



A 



hhh 



1 



1 



1 



2 2/3 + Aiz - /i/i yi + I2h- jJz 2/2 + Aife - /^z 2/iz 

The cross section (j(e+e^ — /i/iZ) is shown for ^/s = 500 GeV at tan/3 = 1.5 as a 
function of the Higgs mass Mh in Fig. 4a. For small masses, the cross section is built up 
almost exclusively hj H ^ hh decays [dashed curve], except close to the point where the 
^Hhh coupling accidentally vanishes (cf. Ref.|P]) and for masses around ~ 90 GeV where 
additional contributions come from the decay A hZ [this range of Mh corresponds 
to Ma values where BR(74 —>■ hZ) is large; cf. Fig.3]. For intermediate masses, the 
resonance contribution is reduced and, in particular above 90 GeV where the decoupling 
limit is approached, the continuum hh production becomes dominant, falling finally down 
to the cross section for double Higgs production in the Standard Model [dashed line]. 
After subtracting the H ^ hh decays [which of course is very difficult], the continuum 
cross section is about 0.5 fb, and is of the same order as the SM. cross section at y/s = 500 
GeV. Very high luminosity is therefore needed to measure the trilinear hhh coupling. At 
higher energies, since the cross section for double Higgs-strahlung scales like 1/s, the 
rates are correspondingly smaller, cf. Fig. 4b. 

Prospects are similar for large tan (3 values. The cascade decay H hh is restricted 
to a small Mh range of less than 70 GeV, with a production cross section of ~ 20 fb at 
a/s = 500 GeV and ~ 3 fb at 1.5 TeV. The continuum cross sections are of the order 
of 0.1 fb at both energies, so that very high luminosities will be needed to measure the 
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continuum cross sections in this case if the background problems can be mastered at all. 

We have repeated the analysis for the continuum process e^e~ —>■ Ahh (cf. Fig.lc). 
However, it turned out that the cross section is built up almost exclusively by resonant 
AH — »• Ahh final states, with a very small continuum contribution, so that the measure- 
ment of the coupling Xhaa is extremely difficult in this process. 

3.2 WlWl hh 

In the effective longitudinal W approximation^ the total cross section for the subprocess 
WlWl — >• hh, Fig. Id, is given by 



+ 




2 \2 

w) 



Hz sin(/3 -a) uz cos(/3 -a) 

^hhh H ; ^Hhh + J- 



1 - yWh 



1 - /i^ 



(16) 



Hz sm{(3 -a) Uz cos(/3 -a) 

^hhh H ; ^Hhh + 



I- Hh 



I- Hh 



[siv?{j3 - a)gi 



+ cos^(/9 - a)g2\ + 



1 



sin^^P — a)g3 + cos^(/? — a)g4 + sin^ 2(/? — a)g5 



with 

91 = 

92 = 

93 = 

94 = 

95 = 



2[iPw - xwPh? + 1 - Pw]lw - ^Phi2Pw - xwPh) 
2ixcf3h - M% + ^PhixcPh - 2(3w) 

Ph[l3hXw{WWw + 14/3^ + 2 - Ww) - ^M'il3lxlr + + 1 - PtvWw + xwVw] 
-[P^ + (1 - + 2/?^ - P^)][lw/xw - yw] - 2Plyw{2Pw - PhXwf 

PhiPhXciSPlxl + UPl,) - APwiSPlxl + P^)][lc + xcvc] 
-Pw[lc/xc - yc] - 2ycPli2Pw - PhXcf 



PhPwl 



w 



x-w X(j 



[2Xwi2Xy^PhPw - XcX^Ph - XcPw) - 2Xy^{Pf^X^ + Py^+l- Py^) 



Xc 



n n ii^Ww + - Pw) + iPlxl + Pwf)] - ^PIMXW + xc) 

PwPh 
PhPwlc 



X. 



c 



X 



[Ax'cPhPw - 2xcXw{Pix'c + Pw + ^-l3w)- ^xUPfix'c + Ptv 



w 



PwPh 



((/?^ + PlxDil - Ptv) + {PWc + Phf)] + WUxcXwPl + 4/?i 



w) 



(17) 



The scaling variables are defined in the same way as before, s^^^ is the cm. energy of the 
subprocess, py^ = {I - AM^/sf/"^ and p^ = {l - AMl/s)^/"^ are the velocities of the W 
and h bosons, and 

xw = {I- 2Hh)/{PwPh) , Xc = (1 - 2Hh + 2hh± - 2hw)/ {PwPh) 

k = \og{xi-l)/{xi + l) , yi = 2/{xl-l) (18) 



*For qualifying comments see footnote 3. 
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The value of the charged Higgs boson mass Mh± in the if^ t-channel exchange diagram 
of Fig. Id is given by Mjj± = M\ + M^. 

In the decouphng hmit, the cross section reduces again to the SAi cross section which 
in terms of gi and g2, defined above, is given by: 



+ 



w) 



hhh 



+ 1 



hhh 



hi 



+ 1 



9i + 
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PwPh 



(19) 



After folding aii with the longitudinal WlWl luminosity [|T^, one obtains the total 
cross section a{e~^e~ v^.v^hh) shown in Fig. 4b as a function of the light Higgs mass 
for tan (3 = 1.5 at ^/s = 1.5 TeV. It is significantly larger than for double Higgs-strahlung 
in the continuum. Again, for very light Higgs masses, most of the events are H —>■ hh 
decays [dashed line]. The continuum hh production is of the same size as pair production 
of SAi Higgs bosons [dotted line] which, as anticipated, is being approached near the 
upper limit of the h mass in the decoupling limit. The size of the continuum hh fusion 
cross section renders this channel more promising than double Higgs-strahlung for the 
measurement of the trilinear hhh coupling. 

For large tan /? values, strong destructive interference effects reduce the cross section 
in the continuum to very small values, of order 10~^ fb, before the SAi cross section is 
reached again in the decoupling limit. As before, the hh final state is almost exclusively 
built up by the resonance H ^ hh decays. 



4. Summa 

It is convenient to summarize our results by presenting Fig.5, which displays the areas of 
the [Ma, tan [3 ] plane in which Xuhh [solid lines, 135° hatching] and \hhh [dashed lines, 45° 
hatching] could eventually be accessible by experiment. The size of these areas is based 
on purely theoretical cuts so that they are expected to shrink if background processes and 
detector effects are taken into account. 

(z) In the case of if — > hh, we require a lower limit of the cross section cr(i7) xBR(if — > 
hh) > 0.5 fb and at the same time for the decay branching ratio 0.1 <BR(if — > hh) < 0.9, 
as discussed earlier. Based on these definitions, Xnhh niay become accessible in two 
disconnected regions denoted by I and II [135° hatched] in Fig.5. For low tan/? , the left 
boundary of Region I is set by LEPl data. The gap between Regions I and II is a result 
of the nearly vanishing Xnth coupling in this strip. The right boundary of Region II is 
due to the overwhelming tt decay mode for heavy H masses, as well as due to the small 
H production cross section. For moderate values of tan (3 , the left boundary of Region I 
is defined by BR(if — * hh) > 0.9. In the area between Regions I and II, H cannot decay 
into two h bosons, i.e. Mh < ^M^. For large tan/? > 10, BR[/7 hh{AA)] is either too 
large or too small, except in a very small strip. Ma — 65 GeV, towards the top of Region 
I. [Note that h and A are nearly mass-degenerate in this area.] 
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(a) The dashed Une in Fig. 5 describes the left boundary of the area [45° hatched] m 
which Xhhh may become accessible; it is defined by the requirement that the continuum 
WlWl — > hh cross section, cTcont, is larger than 0.5 fb. Note that the resonant H ^ hh 
events in Region II must be subtracted in order to extract the Xhhh coupling. 

In conclusion, we have derived the cross sections for the double production of the 
lightest neutral Higgs boson in the J^SSJ^ at e~^e~ coUiders: in the Higgs-strahlung 
process e~^e~ —>■ Zhh, [in the triple Higgs production process e^e~ Ahh], and in the 
WW fusion mechanism. These cross sections are large for resonant H ^ hh decays so 
that the measurement of the triple Higgs coupling Xnhh is expected to be fairly easy for 
H ^ hh decays in the mass range between 150 and 350 GeV for small tan (3 values. 
The continuum processes must be exploited to measure the triple Higgs coupling Xhhh- 
These continuum cross sections, which are of the same size as in the SM., are rather small 
so that high luminosities are needed for the measurement of the triple Higgs coupling Xhhh- 
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Fig.l: Main mechanisms for the double production of the hght MSSM Higgs boson in 
e+e" colhsions: a) e+e" ZH, e+e~ AH and WW — >■ H followed hy H ^ hh; (b) 
e+e- hhZ, (c) e+e" ^ hhA and (d) WW hh. 
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Fig. 2: Cross sections for the production of the heavy CP-even Higgs boson H in e+e" 
coUisions, e'^e~ ZH/AH and e'^e~ — > HvgPg, and for the background process e+e" — > 
Ah [the dashed curve shows | x a^Ah) for clarity of the figures]. The cm. energies are 
chosen = 500 GeV in (a), and 1.5 TeV in (b). 



100 



Fig. 3: The branching ratios of the main decays modes of the heavy CP-even neutral 
Higgs boson H in (a), and of the pseudoscalar Higgs boson A in (b). 



101 



Fig. 4: The cross sections for hh production in the continuum for tan (3 = 1.5: e~^e 
hhZ at a cm. energy of = 500 GeV (a) and WlWl ^ hh at y/s ^ 1.5 TeV (b). 
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Fig. 5: The areas of the [Myi,tan/3 ] plane in which the Higgs self-couphngs \Hhh and 
^hhh could eventually be accessible by experiment at ^/s = 1.5 TeV [see text for further 
discussions] . 
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Loop Induced Higgs Boson Pair Production 
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Institut fiir Theoretische Physik, Universitat Karlsruhe, D-76128 Karlsruhe, Germany. 

Abstract 

We discuss the loop induced production of Higgs boson pairs at high-energy e^e~ 
colliders, both in the Standard Model and in its minimal super symmetric extension. 
The cross sections are rather small, but these processes could be visible with high- 
enough luminosities and if longitudinal polarization is available. 

1. Introduction 

If the genuine supersymmetric particles were too heavy to be kinematically accessible 
in collider experiments, the only way to distinguish between the Standard Model (SM) and 
the hghtest Higgs boson of its minimal extension (MSSM) in the decoupling limit [where 
all the other MSSM Higgs bosons are heavy, and the lightest Higgs boson h has exactly the 
same properties [1] as the SM Higgs boson except that its mass is restricted to be smaller 
than Mh ^ 140 GeV] , is to search for loop induced contributions of the supersymmetric 
particles, which could give rise to sizeable deviations from the predictions of the SM. Well 
known examples of this loop induced processes are the 77 widths of the Higgs particles 
[2] or the process e+e" Z+Higgs which in the MSSM receive extra contributions from 
supersymmetric gaugino and sfermion loops [3]. 

Another type of such discriminating processes is the pair production of Higgs bosons 
which will be analyzed here. In the SM, where it has been first discussed in Ref. [4], 
the process e+e^ H^H^ is mediated only by W and Z boson loops. Fig. la, while 
in the Minimal Supersymmetric extension, additional contributions to the corresponding 
process e^e~ — > hh will originate from chargino, neutralino, selectron and sneutrino loops, 
as well as loops built up by the associated A and bosons; Fig. lb. The cross sections 
for these two processes [as well as for the production of the heavy MSSM Higgs bosons, 
e+e" — > HH^ AA and hH\ have been derived in [5] and here we will summarize the results. 

2. SM Higgs Pair Production 

In the SM, non-zero contributions to the process e^e~ — > H^H^ can only come from 
one-loop diagrams, in the limit of vanishing electron mass. Among these, the diagrams 
involving the one-loop He^e~ vertex [because rrie ~ 0] and those with 7 and Z boson 
s-channel exchanges [because of CP invariance] give zero-contribution; additional contri- 
butions from vertex diagrams involving the quartic WWH^H^/ZZH^H^ couplings are 
proportional rrie and also negligible. The only contribution to Higgs pair production in 
the SM will therefore come from W and Z box diagrams. Fig. la. The expressions of the 
cross sections, allowing for longitudinal polarization of the initial beams are given in [5]. 
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Fig. 1: Feynman diagrams contributing to the Higgs boson pair production process in e'^e 
collisions in the SM (a) and in the MSSM (b). 

The cross sections are shown in Fig. 2 as a function of the Higgs boson mass for two 
center-of-mass energies, ^/s — 500 GeV and 1.5 TeV. Except when approaching the 2Mh 
threshold [and the small dip near the WW threshold], the cross sections are practically 
constant for a given value of the cm. energy, and amount to cr ~ 0.2 fb at ^/s = 500 GeV 
in the unpolarized case. The decrease of the cross sections with increasing center-of-mass 
energy is very mild: at a/s — 1.5 TeV, the cross section is still at the level of cr ~ 0.15 fb 
for Higgs boson masses less than Mh ^ 350 GeV. 

With left-handed polarized electrons, the cross section e^e"*" — > H^H^ is larger by a 
factor of two, while for left-handed electrons and right-handed positrons, the cross section 
e^e^ — s> H^H^ is larger by a factor of four, compared to the unpolarized case. Therefore, 
the availability of longitudinal polarization of the initial beams is very important. With 
integrated luminosities of the order of / £ ~ 100 fb~^ which are expected to be available 
for future high-energy linear colliders, one could expect a few hundred events in the course 
of a few years, if both initial beams can be longitudinally polarized. 

For Mh ^ 140 GeV, the signal will mainly consist of four b quarks in the final state, 
e^e~ H^H^ bbbb, since the dominant decay mode of the Higgs boson in this mass 
range is bb. This calls for very efficient /i-vertex detectors to tag the b jets. Since 

these rare events will be searched for only after the discovery of the Higgs boson in the 
main production processes [5], Mh will be precisely known and the two mass constraints 
m{bb) — Mh, together with the large number of final b quarks, give a reasonable hope 
to experimentally isolate the signals despite of the low rates. For Mh ^ 140 GeV, since 
ffO _^ W^W~ and — > ZZ will be the dominant decay modes of the Higgs boson, 
the signals will consist of four gauge bosons in the final state, e+e^ H^H^ WW, 
leading to eight final fermions. These rather spectacular events should also help to exper- 
imentally isolate the signal. 

3. MSSM Higgs Pair Production 

For the pair production of the light CP-even Higgs boson of the MSSM, e"^e^ — > hh. 
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Fig. 2: The cross sections for, e+e" H^H^, as a function of Mh for ^/s = 500 GeV [dashed 
Hnes] and ^/s = 1.5 TeV [sohd lines]. The lower, middle and upper curves correspond to the 
cross sections with unpolarized, e£ and e^e£ beams respectively. 

several additional diagrams will contribute to the process; Fig. lb. Besides the W and Z 
boson box diagrams, one has the box diagrams with the exchange of the pseudoscalar and 
the charged Higgs bosons, A and and the box diagrams built up by chargino/sneutrino 
and neutralino/selectron loops. The analytical expressions of the cross sections, allowing 
for longitudinal polarization of the initial beams are also given in [5] . 

In Fig. 3, we show the cross section for the process e+e~— as a function of for 
two cm. energies ^/s — 500 GeV and 1.5 TeV and two values of tan /3 = 1.5 and 50. The 
solid lines are for the full cross sections, while the dashed lines are for the cross sections 
without the SUSY contributions. To include the latter we have chosen the parameters 
M2 — —n = 150 GeV, while the common slepton and squark masses are taken to be 
Ml — 300 GeV and Ms = 500 GeV; the parameter At and Af, are set to zero. Only 
the unpolarized cross sections are discussed: as mentioned previously, they are simply 
increased by a factor of 2(4) when the initial beam(s) are longitudinally polarized. 

Let us first discuss the case where the supersymmetric contributions are not included, 
for small tan/3 the cross section is of the same order as the SM cross section and does 
not strongly depend on M^ especially at very high-energies. Although the WWh/Zhh 
couplings are suppressed by sin(/? — a) factors, the suppression is not very strong and the 
W/Z box contributions are not much smaller than in the SM; the diagrams where A/H^ 
are exchanged will give compensating contributions since the hAZ/hH^W couplings are 
proportional to the complementary factor cos(/3 — ct). As in the SM case, the cross sections 
slightly decrease with increasing energy. For large tan values, the factors sin / cos(/9 — «) 
vary widely when M^ is varied. For small M^^ the factor sin(/5 — a) — 0, and the 
contribution of the diagrams with A/H^ exchange dominates. The latter contribution 
decreases with increasing M^ [i.e. with decreasing cos(/3 — a)], until the decoupling limit 
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Fig. 3: The cross sections for e'^e" — hh in the MSSM, as a function of for y/s = 500 
GeV and ^/s = 1.5 TeV and for tan/3 = 1.5 and 50. The soHd curves correspond to the 
full cross sections, while the dashed curves correspond to the cross sections without the SUSY 
contributions. 
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is reached for Mh — 110 GeV. In this case, the factor sm{(3 — a) ^ 1 and the W/Z boson 
loops are not suppressed anymore; one then obtains the SM cross section. 

The contributions of the chargino/selectron and neutrahno/sneutrino loops lead to a 
destructive interference. At high-energies, the supersymmetric boxes practically do not 
contribute; but at low energies, and especially below the decoupling limit, the SUSY 
contributions can be of the order of ~ 10%. We have scanned the SUSY parameter 
space, and the maximum contribution of the SUSY loops that we have found was about 
~ —15%. In the decoupling limit, the SUSY contributions are, at most, of the order of 
a few percent. Because of the rather low production rates, it will therefore be difficult to 
experimentally see this effect. 

4. Summary 

We have discussed the one-loop induced production of Higgs boson pairs at future 
high-energy e~^e~ colliders in the SM and the MSSM. In the SM, the unpolarized cross 
section is rather small, of the order 0.1-0.2 fb. The longitudinal polarization of both the e~ 
and e"*" beams will increase the cross section by a factor of 4. With integrated luminosities 
J jC > 100 fb~^ as expected to be the case for future high-energy hnear colliders, one could 
expect a few hundred events in the course of a few years if longitudinal polarization is 
available. The final states are rather clean, giving a reasonable hope to isolate the signals 
experimentally. In the MSSM, additional contributions to the processes e~^e~ hh come 
from chargino/ neutralino and slepton loops. For hh production, the contributions of the 
supersymmetric loops are in general rather small, being of the order of a few percent; the 
cross sections are therefore of the same order as in the SM. For the processes involving 
heavy Higgs bosons, the cross sections are even smaller than for e"'"e^ hh, and the 
signals will be hard to be detected experimentally. 
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Abstract 



We show that at least one of the Higgs bosons of the Next to Minimal Super- 
symmetric Standard Model can be detected at future Linear Colliders of 500, 1000 
and 2000 GeV cm. energies. 

1. Introduction 

The Next to Minimal Supersymmetric Standard Model (NMSSM) [1-3] is a minimal 
extension to the Minimal Supersymmetric Standard Model The NMSSM provides 
the most economic solution of the so called /i-problem of the MSSM by introducing an 
additional Higgs singlet superfield M = {N, ip]\f, F^r) with a Higgs singlet A^, a higgsino 
singlet ipN and an auxiliary field F/y. Together with the two Higgs doublets superfields 
TCi^2 = (-f^i,2, ^^1,2, -^1,2) the superpotential of the NMSSM is given by 



where Ax and A^. are soft breaking mass parameters. 

Hi,H2 and develop vacuum expectation values Vi,V2 and x respectively. The 
NMSSM contains three scalar Higgs bosons 5*1, S2 and S3 with masses ms^ < ms^ < ms,^, 
two pseudoscalar Higgs bosons Pi and P2 with masses mp^ < mp^ and a charged one 
with mass me- 

The Higgs sector has 6 free parameters A, k, tan j3 = f 1/^2, x, Ax and A^. 
2. Constraints on parameters 

A remarkable result of the MSSM is the tree level bound of < mz cos 2/3. This is due 
to the fact that all quartic terms have gauge coupling constants. In case of the NMSSM 
there is a quartic term with the coupling constant A. It turns out that the upper bound 
of A may be relevant for m^^ . An effective way of determining this bound is RG-analysis 




(1) 



The soft breaking part of the Higgs sector is given by 




(2) 



[1-8]. 
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The one loop RG-equation of A is coupled with that of k and ht, the Yukawa coupling 
constant of the top quark (neglecting other quarks) and is given by 

where t = In /i and /i being the renormalization scale. By demanding no Landau pole up 
to the GUT scale one can determine from eq. (|) the upper bound of A and k and the 
lower bound of tan/5 at the electroweak scale. We plot our results in Fig. 1 (Fig. 2) for 
rrit = 175 GeV (190 GeV). They show that Xmax decreases with increasing k^. The lower 
bound of tan (3 is about 1.24 for rrit = 175GeV and 2.6 for rrit = 190GeV. For tan/? > 3 
Xmax is almost independent on tan/?. The upper bound of k is about 0.7. 

3. Mass upper bounds 

The tree level bound of is given by 

f 2A^ \ 

< m% cos^ 23 H — ^ k sin^ 23 = mimax (6) 

" V 9i+ 92 '^y ^1 ^ ^ 

The upper bound of ms^ and m^g can be expressed in terms of 1713^^^^ and ms^ 



2 . 2 _ "^1— - {Rl + RlWs, 



1 - (i?? + 



where i?j = f/ji cos j3 + ?7j2 sin (3 and ?7jj is the 3x3 orthogonal matrix which diagonalizes 
the scalar mass matrix. 

Rl and R2 satisfy the unitarity condition < i?^ + r| < 1. The tree level upper bound 



(D yields ms, < mz for A^ < {gf + gi)/2 = (0.52)^ and ms, < pl{gl + gl)Xmz = 
1.92Xmz for A^ > (0.52)^. Using Xmax = 0.64 — 0.74 from section 2 the tree level relation 
yields mg, < 113GeV - 131GeV. 

As in the case of the MSSM the contributions of radiative corrections may change 
this result considerably. Several groups calculated higher order contributions to the mass 
matrices using the one loop effective potential and determined the corrected upper bound 



*We obtain Xmax = 0.64 - 0.74 for m* = 175 - 190 GeV. 
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The result in our notation |Q] is given by 



^li < rn\ [ cos^ 2(3 + — 2 ^^^^ 2(3] + a cos^ /3 + /3 sin 2/5 + 7 sin^ /3 



+ 92 



(9) 



with (At = -At + Ax cot /?) 



a 



7 



1 / XxAt V / ^ A ^ 



2 

/ m," \ 

-XxAt 
ill KmiVi 



y ^ , aat 



6m| 



57r \ V^ 



2 In 



mi 



rrir 



6m^ 



(10) 
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(12) 



In this result only top and stop contributions were taken into account. We numerical 
calculated m^mai, in the region 175GeV < rrit < 190GeV, 250GeV < x, Ax, At, mi < 
lOOOGeV and 2 < tan/5 < 20 and obtained || 



120 < ms^-- < 156GeV 



(13) 



4. Production of scalar Higgs bosons at 



e ' e 



Colliders 



The upper bound ms^ < 120 — 156GeV suggests that the accessible area of the parameter 
space at LEPl with ^/s = rriz might be very small. Actually we showed that the existing 
LEPl data do not exclude the existence of Si with m^^ = OGeV |]12|. 

For colliders with ^/s = 500, 1000 or 2000GeV the situation is different. In this case 
the production cross section of one Si via the Higgsstrahlung e~^e~ —>■ ZSi with real Z 
and 5", is always possible as the collider energy is larger than Et = 212 — 248GeV. Et is 
a kind of threshold energy and is an important quantity of a model. 

In this case it is possible to derive a lower bound for the production cross section cxj 
of Si as a function of the collision energy only. This lower bound would give information 
about how far the model could be tested. 

In order to derive the lower bound of ai we consider the production cross sections of 
^i, 5*2, Ss via the Higgsstrahlungs process, denoted by ai, (T2, as, which can be expressed 
in terms of the standard model Higgs production cross section asM and -Ri and R2 defined 
in section 2: 



0-1 (ms^; 

A useful observation is that ai(m 



RlcTsMifUs,) 

RlasAiims^) 

(1 - ^1 - Rl)(^SMims:^ 



(14) 
(15) 
(16) 



5, 



< ai{msi) which allows one to derive parameter 



independent lower limits on 0", as we will see in the following. 



Ill 



First we determine at a fixed set of ms^^, Ri, R2 the cross sections ai{Ri, R2,msi), 
(72(i?i, i?2, ^s™^) and (T3(i?i, i?2, ^5™^)- Then we keep Ri and R2 fixed, but vary ms^ 
from its minimum to its maximum value and determine the quantity a{Ri,R2) defined 

by 

cr(i?i,i?2)= min [max(cri, cr2, cxs)] (17) 

where ai = ai{Ri, R2,ms^), a2 = a2{Ri-, R2-,rns^'^^) and cxa = a3{Ri, R2,ms^'^^). As a 
last step we vary R^ and i?| from to 1 with i?^ + < 1 and plot a{Ri, R2) in the Ri- 
i?2-plane. It is plausible that cr(i?i, R2) = for ^/s < Et = + ms-^^^ = 212 — 248GeV 
in the entire i?i-i?2-plane. This is the case for LEP2 with a/s < 205 GeV. Therefore this 
method docs not give any results for LEP2. For ^/s > Et cr{Ri, R2) never vanishes and 
the minimum of cr(i?i, R2) is a parameter independent lower limit of one of the a^. This 
minimum is thus a characteristic quantity of the model. 

In Fig. 3 we plotted a{Ri,R2) for = 500 GeV and rngma. = 145 GeV. The 
minimum is about 16 fb. When the discovery limit is about 30 events, one would need 
a luminosity of about 25 fb, which is a realistic one. Fig. 4 (Fig. 5) shows o"(-Ri,i?2) 
for ^/s = 1000(2000) GeV with minimum cross section of 4 fb (1 fb). Fig. 6 shows 
the minimum of a{Ri,R2) as a function of y/s and ms^<^^ as a parameter. We see that 
the effect of msmax on amin is very big around ^/s — 300 GeV, but rather small for 
v/s> 500 GeV.' 

Fig. 7 shows the tree level cross sections (T2 and for an exemplary set of param- 
eters with the contributions from (i) the Higgsstrahlungsprocess e+e^ —>■ ZSi —>■ bbSi, (ii) 
the process where Si is radiated off from 6 or and (ni) e+e" — > Z — > PjSi — > bbSi, where 
Pj{j — 1, 2) is a pseudoscalar Higgs boson. 

Fig. 8 shows the same as Fig. 7, but with one loop contributions via the effective 
potential. The higher order contribution is rather important for the energy region around 
150 GeV and decreases with ^/s. In this parameter region the dominant production is 
that of 5*2 at y/s = 500 GeV and is about 13 fb. 

We conclude that the Higgs sector of the NMSSM can most probably be tested con- 
clusively at the future hnear e''"e~-colliders with 500, 1000 or 2000 GeV cm. energies. 
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Figure 1: Upper bound for X as a funtion of tan (3 for rritop = 175 GeV. 



Figure 2: Upper hound for X as a funtion of tan /3 for rritop — 190 GeV. 
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Figure 3: a{Ri, R2) as definded in the text for y/s — 500 GeV. 



Figure 4: a{Ri, R2) for y/s = 1000 GeV. 
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Figure 5: a{Ri, R2) for y/s = 2000 GeV. 



Figure 6: Minimal value of a{Ri, R2) as a function of \/s for various values of ms^""^ 



Figure 7: Cross section for e^e~ Zbb for A\ = 220 GeV,Ak = 160 GeV,x = 
1000 GeF, tan/3 = 2,k = 0.04, A = 0.12. Masses and mixing angles have calculated 
from the tree level protential. 



Figure 8: The same as above, but with masses and mixing angles obtained from the one 
loop effective potential. The top mass is 1 75 Ge V, mx — mi — 1 TeV and At — Ge V. 
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